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ABSTRACT 


The  Helmholtz  integral  and  the  Kirchhoff  approximation  have  been 
used  to  develop  predictions  for  the  transmission  of  sound  through  a 
rough  air-sea  interface.  A  model  study  was  conducted  with  wind-driven 
surfaces  generated  in  a  large  anechoic  tank.  Root  mean  square  wave 
heights,  a,  ranged  from  .05  to  .41  cm,  windward  correlation  lengths,  L  , 
from  0.8  to  4.2  cm.  The  frequencies  used  (10-43  kHz)  were  scaled  to 
be  equivalent  to  low  audio  frequencies  with  moderate  seas.  Although 
the  coherent  component  of  the  transmitted  acoustic  intensity  showed  an 

exponential  decrease  with  increasing  values  of  the  roughness  parameter, 

2  2  2 

R  [  =  kp  a     (c2/c-j  cos  e,  -  cos  e2)  ],  the  decrease  was  greater  than 

predicted.  (Subscript  2  refers  to  propagation  constant,  speed  and 
angle  of  transmission  in  water;  subscript  1,  in  air.)  This  lack  of 
agreement  appeared  to  be  caused  by  violation  of  the  Kirchhoff  require- 
ment, L/x2  >>  1 .  An  empirical  correction  factor  a  =  1  +  4.8e  u  ' 
was  determined  to  give  a  corrected  roughness,  aR.  The  corrected  theory 
appears  to  be  valid  for  both  coherent  and  incoherent  components  of 
intensity  for  all  values  of  L/x  provided  that  aR  <  4.0  .  For  higher 
values  of  roughness  there  is  a  "leveling  off"  of  the  transmission  loss 
of  the  coherent  component  similar  to  that  recently  observed  but  unexplained 
in  rough  surface  scattering. 
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I.   INTRODUCTION 

The  problem  of  transmission  of  airborne  sound  into  the  ocean  has 
received  very   little  attention.  Only  recently  has  it  become  apparent 
that  a  more  detailed  study  of  the  problem  than  has  been  made  in  the  past 
should  be  undertaken.  Most  of  the  theoretical  studies  that  have  been 
done  in  this  area  have  dealt  with  a  smooth  air-sea  interface  and  the 
rough  air-water  transmission  studies  reported  to  date  have  been  with 
deterministic  surfaces  [1],  The  work  reported  in  this  thesis  is  part  of 
a  project  in  progress  at  the  Naval  Postgraduate  School  that  has  included 
a  series  of  ocean-acoustic  experiments  as  well  as  the  theoretical  and 
model  studies. 

The  author  has  had  the  satisfaction  of  participating  in  two  success- 
ful ocean  experiments  with  Professor  H.  Medwin  and  LCDR  R.  A.  He! big. 
LCDR  He! big  will  report  on  this  phase  of  the  project  when  the  analysis 
has  been  completed.  These  experiments  were  conducted  from  the  research 
platform  FLIP  operated  by  the  Scripps  Oceanographic  Institution,  and 
involved  the  recording  of  the  underwater  sound  fields  of  several  opera- 
tional Navy  aircraft.  Ocean  wave  data  were  collected  simultaneously 
with  the  acoustic  data  collection  as  were  the  other  important  environ- 
mental parameters.  It  is  hoped  that  the  results  of  these  experiments  can 
be  better  interpreted  on  the  basis  of  the  theoretical  and  model  sea 
studies  reported  here. 

The  objective  of  this  thesis  has  been  to  reduce  the  sound  transmis- 
sion problem  to  the  simplest  useful  case  and  to  apply  the  Helmholtz 
integral  in  the  Kirchhoff  approximation  to  develop  a  theory  that  parallels 
that  derived  for  rough  surface  scattering.  The  rough  surface  scattering 


problem  has  received  a  great  deal  of  attention  and  this  method  of  approach 
has  given  reasonably  good  results  for  the  scattering  of  underwater  sound 
from  a  realistically  rough  sea  surface.  Apparently  no  one  has  tried  to 
apply  this  technique  to  a  rough  surface  transmission  problem  although 
there  seems  to  be  no  reason  why  this  can  not  be  done. 


II.     THEORY 

The  development  given  here  applies  the  Helmholtz  integral   and  the 
Kirchhoff  approximation  to  describe  the  transmission  of  acoustic  waves 
through  a  randomly  rough  air-sea  interface.     The  formulation  of  this 
problem  parallels   that  used  for  acoustic  scattering  from  a  rough 
surface  by  Medwin  [2].     The  monograph  of  Beckmann  and  Spizzichino  [3] 
presents  a  detailed  discussion  of  the  application  of  the  Helmholtz 
integral   to  the  scattering  of  electromagnetic  waves  that,  with  minor 
modifications,  is  directly  applicable  to  the  problem  of  acoustic  trans- 
mission through  a  rough  interface. 

A.     AN   INTEGRAL  EXPRESSION   FOR  THE   RECEIVED  ACOUSTIC  PRESSURE 
The  Helmholtz  integral, 


p2  ■  £  s'<p-& -*!&>*  o) 


expresses  the  complex  acoustic  field,  p«,  at  a  receiver  in  the  far  field 

of  an  ensonified  area,  S,  in  terms  of  the  field  and  its  normal  derivative 

ikr 
on  S  and  the  scalar  "Green's  Function",  ¥  =  — —  .  A  single  frequency 

harmonic  wave  is  assumed  and  the  time  dependence  suppressed.  For  conven- 
ience the  quantity  p2  is  described  as  the  "received  acoustic  pressure" 
with  the  understanding  that  it  represents  the  amplitude  and  phase  of  the 
received  acoustic  wave.  The  exact  value  of  the  field  and  its  normal 
derivative  on  S  is  in  general  unknown  and  the  Kirchhoff  method  consists 
of  approximating  these  values  by  those  that  would  exist  on  a  tangent 
plane  at  that  point.  Within  this  approximation  the  boundary  conditions 
can  be  stated  in  terms  of  the  incident  field,  p,  , 


p     =     p1  +  Rp1     =     Tp1 


9p,  3Rp. 

W    '     ~3rT 


^    =  =    i(l  -  RKt-iOPi  '(2) 


where  R  and  T  are  the  reflection  and  transmission  coefficients  defined 
for  infinite  planes. 

.    Before  proceeding  with  the  evaluation  of  (1)   it  is  convenient  at 
this  point  to  specialize  and  introduce  the  reflection  coefficient  appro- 
priate for  acoustic  transmission  from  air  to  water.     This  results  in  a 
simplification  of  the  integral   since  R  may  be  set  equal   to  one  and  the 
normal   derivative  of  the  field  vanishes  at  the  boundary.     It  is  at  this 
point  that  the  development  first  diverges  from  that  for  underwater 
scatter,  where     R  =  -1     and  the  field  itself  must  vanish  at  the  boundary, 
The  Helmholtz  integral   can  now  be  written  in  a  simplified  form. 

p2     =     h    f     2  pl    (-ik2-nH0   dS  (3) 

Figure  1  is  a  sketch  of  the  geometry  of  the  problem  that  shows  the 
coordinates  that  will  be  used  and  labels  the  angles  that  are  defined. 
A  list  of  the  notation  to  be  used  follows. 

p1A  =  acoustic  pressure  at  unit  distance  from  source 

p1n  =  -~  eUK"|r-|-u)t)^  incident  acoustic  pressure  at  origin 

i  (t     r) 
p,   =  p,Q  e  v  1  '        ,  incident  acoustic  pressure  at  element  dS 

P2   =  acoustic  pressure  at  the  receiver 

dS   =  dxdy/cos  y   ,  the  area  of  a  surface  element 

C    =  C(x,y)  ,  the  surface  height  at  (x,y) 

n    =  i,(-sin  ycos  a)  +  i2(-sin  y  sin  a)  +  i3  cos  y   ,  normal  to 
surface  element 
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see  below 


FIGURE  1.  THE  GEOMETRY  OF  THE  PROBLEM 
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Y        =     angle  the  normal  makes  with  the  z  axis 

a        =     angle  the  projection  of  the  normal   onto  the  xy-plane  makes 
with  the  x  axis 

^1       =     -£-  (i-|   sin  e-|  +  ig  cos  6-j)    ,  vector  propagation  constant 

of  incident  field 

lL  =  -jr-  (i-i   sin  92  cos  9^  +  io  s"i'n  eo  sin  e3  +  ^3  cos  e2^> 

vector  propagation  constant  of  received  field 

9,  =  angle  of  incidence 

&2  =  angle  that     k2     makes  with  the  z-axis 

03  ■  angle  that     lL     makes  with  the  xz-plane 

r  =  i-,x  +  i^y  +  i"3^   ,  the  radius  vector  to  the  surface  element 

r,  =  distance  from  source  to  origin  of  ensonified  area 

r, '  =  distance  from  source  to  surface  element 


r?       =     distance  from  origin  of  ensonified  area  to  receiver 

e  element  to  receiver 

e"1   2'r       t  Green's   Function 


r2'     =     distance  from  surface  element  to  receiver 
eik2r2'  e-k2r2       _tf    + 


r2  r2 


(I) 


=     -  ik9-ny 
S  d 


Continuing,  equation  (3)   can  be  rewritten, 


i2p10eik2r2 


p        =      , f    k 

C  47rr2  S 


/t,.neiK'r     dS  (4) 


where 


K  =   ki  "  4  =   h\  +  12^  +  hh 


K      =     k,sin  0,   -  k2sin  62  cos  93 
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and 


K   =  -kp  sin  0p  sin  03 


K   =  k,  cos  e,  -  k«  cos  9p 


]L-n  =  kp  cos  62  [-tan  6^  sin  y   cos  a  cos  93 


tan 


»2  sin  y  sin  a  sin  9^  +  cos  y] 


Assuming  uniform  ensonification  over  a  rectangular  area  2X  by  2Y, 


-i2pin  elk2r2  k9   cos  Q9  Y   X 

P2     - /      /     f(elte?fevY,a)e,N"rdx4y     (5) 

c  4Trr„  -Y     -X  '     L     J 


.it? 


where 


f     =     -tan  92  tan  y  cos  a  cos  e~  -  tan  6p  tan  y  sin  a  sin  63  +  1    . 


The  smooth  surface  transmitted  pressure  has  a  maximum  value  p?Q,  where 
-i2p10e^2 


p20     = 


A  cos 


4irr, 


A  =  4XY 


(6) 


The  integral   of  (5)  may  be  written  using  (6), 


^20 
A 


./  ./  < a  §  +  b  %  + 1  >  ^r  dxdy 


(7) 


where  ■£    =  tan  y   cos  a 


■^    =  tan  y  sin  a 
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and  a    =     -tan  e2  cos  e3 


b    =     -tan  e2  sin  63 


Now  there  are  three  integrals  to  evaluate  two  of  which  must  be  integrated 
by  parts.     Starting  with  the  first  term 


//    iff    e1^  dx  dy 


(8) 


let  u     =     eiKxx 


du     =     iK^  eiK*x  dx 


dv    =     iKzeiKzc    |f-      dx       ,  v     =    eiKz^ 


Then  since 


/  udv    =     uv  -  /  vdu 


(8)  can  be  evaluated 


•$-    /  eV  /  eiKxx     (iKzeiKz^     f§    dx)       dy 


iK. 


/  e1*  V  E  eiKxX  eiKz^ 


/  eiKzc  eiKxx  iKY  dx  ]  dy 
-X  -X 


iK 


a       [      /  eiKxx 


z         -Y 


_x    M*    e'Vdy-iVy/_x/e'^dxdy] 


Similarly  for  the  second  term 


%    f* 


iK-r 


Y  bKv         Y       X       lit-? 

Y  dx     -     -J-        f      f     elK  r       dxdy 

Kz     -Y     -X 


(9) 
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combining  the  three  integrals  gives, 


F        !       *     Jfrr 


P2     =     p20  S"    _/  _/     e  dxdy  +  B  (10) 


■Y     -X 

where 

aK     +  bK 
F    =    f(en,e9,eJ    =    1  -  ( 


'1>«2,03'  '        v         Kz 


F(ere2,e3)    =    l 


(c2/c-|   sin  0-|   -  sin  9„)   tan  62  cos  93 
c2/c,   cos  9,   -  cos  92 


2 

sin  62  tan  92  sin  e3 


c2/c-|  cos  0,  -  cos  e 


It  can  be  seen  that  F  reduces  to  unity  when  the  receiver  is  in  the 
direction  predicted  by  Snell's  law  of  refraction,  that  is  when 
c2/c,  sin  9,  -  sin  e?  =  0  and  e3  =  0.  In  addition,  F  diverges  signifi- 
cantly from  unity  only  for  angles  far  from  the  Snell  angle  for  given  angle 
of  incidence.  The  behavior  of  F(9-,,92,0)  is  given  in  Table  1,  which  may 
be  compared  with  Table  3.1  in  Ref.  [3]  that  gives  values  for  a  similar 
function  encountered  in  the  surface  scattering  problem.  A  difference  that 
is  observed  is  the  appearance  of  negative  values  of  F  for  refracted  angles 
far  from  the  Snell  angle  and  nearly  horizontal.  The  acoustic  intensity 
in  these  directions  will  be  relatively  low,  since  other  factors  will 
dominate  the  behavior  of  the  transmitted  field,  and  it  will  be  difficult 
to  determine  experimentally  if  this  effect  apparently  introduced  by  F  is, 
in  fact,  real . 
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TABLE  1 
Values  of  F(e, ,e?,0) 


F(e1,e2,0) 


(c2/c.j  sin  e,  -  sin  e2)  tan  e2 


Cp/c,  cos  e,  -  cos  e„ 


^^O 

0° 

3° 

6° 

9° 

12° 

-80 

-0.30 

-0.61 

-0.93 

-0.13 

-1.59 

-70 

0.37 

0.21 

0.55 

-o.n 

-0.27 

-60 

0.62 

0.52 

0.41 

0.31 

0.19 

-50 

0.76 

0.68 

0.61 

0.54 

0.46 

-40 

0.85 

0.80 

0.75 

0.69 

0.63 

-30 

0.92 

0.88 

0.84 

0.80 

0.76 

-20 

0.96 

0.94 

0.92 

0.89 

0.86 

-10 

0.99 

0.98 

0.97 

0.96 

0.94 

0 

1.00 

1.00 

1.00 

1.00 

1.00 

10 

0.99 

1.00 

1.02 

1.03 

1.04 

20 

0.96 

0.99 

1.01 

1.04 

1.06 

30 

0.92 

0.96 

0.99 

1.03 

1.07 

40 

0.85 

0.91 

0.96 

1.01 

1.07 

50 

0.76 

0.83 

0.91 

0.98 

1.05 

60 

0.62 

0.72 

0.82 

0.92 

1.03 

70 

0.37 

0.53 

0.68 

0.84 

0.99 

80 

-0.30 

0.28 

0.31 

0.61 

0.92 

16 


The  second  term  in  (10)  represents  an  "edge  effect"  and  can  be  shown 
to  be  small  compared  to  the  remaining  integral  if  the  ensonified  area  is 
many  wave  lengths  in  dimension. 


»   -   P20   r   3 


Y  it.-? 

/  elK  r 


-X       lKz  -X 


-Y 


dx  ]    (11) 


To  see  this,  it  is  argued  that  because  the  conditions  at  the  periphery 
of  a  large  ensonified  area  can  not  be  a  major  factor  in  the  overall  result 
C  can  be  set  equal  to  zero  at  the  edges  without  significantly  altering 
the  predicted  field.  With  this  assumption  the  integration  can  be  per- 
formed and  gives 


aK  +  bK 

-X-   )  sine  KxX  sine  k  Y 


=  p20  (  ""¥ 


or 


B  =  p2Q  (1  -  F)  sine  KxX  sine  k  Y 


02) 


It  can  be  observed  that  because  of  the  sine  terms,  B  can  only  have  a 
significant  value  within  a  narrow  wedge  centered  on  Snell's  direction. 
But  in  this  direction  F  is  approximately  equal  to  unity;  so  we  have 
B  «  p20  .  On  the  basis  of  the  above  argument  the  edge  contribution  is 
dropped  and  (10)  becomes 


p20 


Y   X 
/   / 
-Y  -X 


e1*'*  dxdy 


(13) 


Equation  (13)  is  an  integral  expression  for  the  received  acoustic 
pressure.  The  form  of  this  equation  is  identical  to  that  developed  for 
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the  surface  scattering  problem.  The  steps  from  this  point  on  are  iden- 
tical to  those  in  surface  scattering  but,  of  course,  the  results  will 
be  somewhat  different  since  the  terms  in  (13)  reflect  the  different 
geometry  and  refraction  encountered  in  rough  surface  transmission.  To 
procede  a  further  specialization  will  be  made  by  specifying  a  randomly 
rough  surface  and  the  first  and  second  moments  of  the  transmitted 
acoustic  signal  will  be  derived. 

B.  ASSUMPTIONS  AND  APPROXIMATIONS 

In  order  to  derive  (13)  it  has  been  necessary  to  make  a  number  of 
assumptions  and  approximations  which  may  limit  the  applicability  of  the 
theory.  To  summarize,  it  has  been  assumed: 

1)  The  incident  wave  is  plane  and  ensonifies  a  limited  area. 

2)  The  receiver  is  sufficiently  far  from  the  surface  to  regard  the 
received  waves  as  plane,  i.e.,  the  "far  field  approximation"  is 
made. 

3)  Mutual  interaction  of  the  surface  irregularities  may  be  neglected. 

4)  The  ensonified  area  is  large  compared  to  an  acoustic  wavelength. 

5)  The  surface  is  considered  perfectly  rigid  such  that  the  reflection 
coefficient,  R  =  1 . 

Qyfi    The  Kirchhoff  approximation  is  valid. 

Of  the  above,  the  last  is  probably  the  most  serious  limitation  to 
the  validity  of  the  solution.  The  Kirchhoff  approximation  has  been 
noted  by  Beckmann  and  others  to  require  that  the  minimum  radius  of  cur- 
vature of  the  surface  be  large  compared  with  the  wavelength  of  the 
incident  radiation.  For  the  transmission  problem  there  are  two  wave- 
lengths at  the  interface  and  it  seems  reasonable  that  the  reradiation 
wavelength  (x2),  which  in  this  case  is  the  larger  wavelength,  must  be 
chosen  as  the  critical  one.  The  requirement  of  the  Kirchhoff  approxima- 
tion seems  to  limit  the  application  of  (13)  to  a  certain  class  of 
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surfaces.  The  surface  must  be  a  gentle  one  without  a  large  amount  of 
sharp  edges  or  irregularities  with  small  radii  of  curvature.  For  a 
realistic  wind-driven  sea  this  requirement  can  be  expressed  in  terms  of 
a  characteristic  length  that  may  be  derived  from  the  two-dimensional 
spatial  correlation  or  equivalently  from  the  sea  surface  directional 
spectra. 

C."  THE  MEAN  ACOUSTIC  PRESSURE 

The  average  received  acoustic  pressure  can  be  calculated  for  a  sta- 
tistically rough  surface  by  taking  the  space  average  of  equation  (13). 


<P2>=  P20£   ./-/  *W  +  K/)<^>  ^  ™ 


It  has  been  shown  that  the  wave  height  distribution  of  ocean  waves  is 
very  nearly  Gaussian.     Using  a  Gaussian  probability  density  function 

w(?)     =    _J^      e"A2°2 

cW2tt 

7 
Equation  (^f  can  now  be  evaluated. 

<(p2)    =  p20Fe~R/2   sinc  KxX  sinc  KyY  (15) 


where    R  =  Kz2a2 


The  significant  parameter  in  equation  (7)  is  a  "roughness  parameter' 

similar  to  one  encountered  in  rough  surface  scattering  and  denoted  by 

2  2 

Beckmann  as  g  =  K  a  .  Comparison  of  R  and  g  shows  some  interesting 

differences: 


19 


2  2  2 

R  =  k2  a     (c2/c-j  cos  e1  -  COS  6  ) 

g  =  k/o2  (cos  Q]   +  cos  e2)2 

A  significant  difference  is  that  the  strong  dependence  on  geometry  seen 
in  g  is  almost  completely  absent  in  R.  While  g  can  vanish  as  the  angles 
approach  grazing,  the  geometry  dependence  in  R  is  dominated  by  the  term, 
c2/c,  cos  e,  and  this  quantity  must  remain  nearly  constant  because  of 
the  limited  range  of  incident  angles,  e,  ,  that  will  result  in  significant 
sound  transmission  for  sound  traveling  from  air  into  water.  If  R  is 
evaluated  for  those  refracted  rays  predicted  by  Snell's  Law,  it  can  be 
seen  to  be  a  minimum  at  normal  incidence  and  to  increase  gradually  as 
the  angle  of  incidence  increases.   It  is  also  apparent  that,  for  the  same 
frequency  and  mean  square  wave  height,  the  surface  appears  "rougher"  for 
sound  transmission  from  air  into  the  sea  than  it  does  forthe  scattering 
of  underwater  sound.  R  will  be  a  minimum  of  three  times  greater  than 
g  at  normal  incidence  and,  of  course,  much  larger  if  the  practically 
important  case  of  near  grazing  incidence  for  surface  scatter  is  compared 
to  nearly  horizontal  transmission. 

D.  THE  VARIANCE  OF  THE  ACOUSTIC  PRESSURE 

In  order  to  evaluate  the  mean  received  intensity,  the  variance  of  the 
acoustic  pressure  must  be  determined  in  addition  to  the  mean  acoustic 
pressure.  The  variance  of  p2  can  be  evaluated  through  the  formula 

V<P2)  =  <p2p2*>  -  (P2><P^ 

which  can  be  expressed 
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V(Po)  ■   P, 


2  ,  F  x2 


(  I)2  ffff  euKx(x-x')  +  Ky(y-y)) 


'T         K20   v  A 

[<(e'KzU-l'Y)>     .    <(eiKz^>  ]  dxdydx'dy'         (16) 

For  a  Gaussian  distribution  of  heights, 

(VMc-c'))    =   e-R  +  RC 

where  C  is  the  correlation  coefficient. 

It  is  now  assumed  that  the  surface  is  isotropically  rough  and  a  Gaussian 
correlation  function  is  specified.  Although  this  assumption  is  quite 
obviously  incorrect  for  a  realistic  sea,  it  has  yielded  some  useful 
results  in  the  surface  scattering  problem  and  it  simplifies  the  mathe- 
matics considerably. 

The  sea  surface  spatial  correlation  is  assumed  to  be 

C(0  =  e"£  /L  (17) 

Transforming  (16)  to  polar  coordinates,  using  the  assumption  L  <<  A. 

V{p2}  =  p2()2£:  ;"     )     e1  (Kxcos*+Kysin*'  [e-R(1-C)-e-R]  d*£d.  (18) 
Using  the  series  expansion  for  the  exponential  the  bracketed  term  becomes 


[e-«  \      t.    e-"*2/L2] 


Integrating  over  i>   gives 
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V<p2)  =  2i  p2„2  £-   /"  J0(«Kxy)  [e"R  z   ^e-ni/L]  ldi  (19) 


where  Kxy  =   K,2  ♦  Ky2 


Uniform  convergence  permits  the  change  of  order  of  summation  and  inte- 
gration. The  resulting  form 


v{p2,  =  2,  p2Q2  F_  e-R  j^  £   /oo(Kx/)  iM2/L2   id(  (20) 
is  a  standard  integration  which  gives 

ton  )     n   2F2  ^   P"R   V   R"   /V2^4"  ,„, 

V{p2)  -  P20  F   T  e    ^  J,,,-  e  (21) 


It  is  possible  to  get  less  complicated  expressions  than  (21)  for 

three  cases:  a)  R  <<  1   b)  R  >>  1   c)  K   =  0  ,  where  K   =0 

xy  xy 

corresponds  to  transmission  in  Snell's  direction. 

a)     For  R  <<  1  we  may  take  only  the  first  term  of  the  series  and  obtain 

2     2     ttL2  -R     "Kxv2l2/4 

V(P2)   -   P2q     F^ReRe     xy  (22) 


b)     For  R  >>  1   the  mean  acoustic  pressure  approaches  zero  and 
V<P2*  0D2P2*/'  '     Ecluation   08)  becomes 

p     2 
<P2D2*>    -     2,     -J*-    F2     /W)  e"R(1   "  C)       ^ 
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or  using  the  Gaussian  correlation  coefficient 

<y2.> .  2,  ?4- p2  /j0(v)e"RCi"exp("2/L2)]  <d*  <23> 


It  can  be  shown  that  for  R  >>  1  this  integral  receives  significant 
contributions  only  from  the  neighborhood  of  i   =  0.  Therefore,  using 
the  first  two  terms  of  the  expansion  of  the  correlation  function, 
CU)     1  -  U2/L2)  ,  (23)  becomes 

^p2p2*>   =   2tt  -|9_  F2   /  0o(Kxy£)  exp  (-FU2/L2)   Ad*      (24) 

and  we  can  evaluate  this  integral 


<P2P2*>      =     P202p' 


exp(-K  „  LV4R)  (25) 


AR   CA^  "xy 


Using  the  expression  for  the  mean  square  slope  of  the  rough  surface, 
i      =  -2y  ,  (25)  can  be  written,  again  for  R  >>  1 


2  K 

V{p2>~  <^P2P2*\    -   P202  ^o—o  exp  (  -    y2  )       (26) 
£  V  *■/  R»l     ^U   K^Ar         2K  V 


c)  A  special  case  of  interest  is  an  evaluation  of  the  variance  of  the 

acoustic  pressure  in  Snell's  direction,  where  K   =  0.  Using  equation  (21) 

xy 


v<p2>  .  p202  4-  ( e"R  j,    £r )  <27> 
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Medwin  [2]  has  evaluated  the  function,  S(R),  where 


S(R)  =  e  R  I      ^r  (28) 

n=l 


A  plot  of  S(R)  is  given  in  Ref.  [2].  It  can  be  noted  that  for  R  <<  1 , 

S(R)  ■*-  R  and  for  R  >>  1  ,  S(R)  ■+  1/R  .  Using  these  asymptotic  values 

and  K   =  0  ,  equation  (27)  agrees  with  the  low  roughness  case  stated 
xy 

in  (22)  and  with  the  high  roughness  case  given  in  (25). 

E.  THE  MEAN  SQUARE  ACOUSTIC  PRESSURE 

It  is  now  possible  to  write  the  general  expression  for  the  mean 
squared  acoustic  pressure. 


<P2P2*)>   =  P202  p2  e"R  [  (sinc  KxX  sinc  KyY)2  + 


n=l 
for  low  roughness,  R  <<  1  ,  we  have 

2  -2  -R 


4-    1       HnV  exp  (-Kxy2L2/4n)  ]    (29) 


<(P2P2*>   =  p20  F  e   [  (sinc  KxX  sinc  KyY)   + 


for  high  roughness,  R  >>  1  ,  we  have 


<e2P2*>  ■  p2o2  f2  ^7  «P  <  ^2  > 
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R  exp(-K  L2/4)  ]  (30) 


(31) 


The  first  term  in  (29)  represents  the  coherent  component  of  the  trans- 
mitted signal,  while  the  second  term  is  the  incoherent  component.  Since 

p 

L  «   A  ,  the  coherent  component  dominates  for  R  <  1.  But  as  R  increases 

the  coherent  term  fades  rapidly  until  we  have,  for  R  >>  1 ,  a  completely 
incoherent  signal . 

F.  STATISTICAL  DISTRIBUTION  OF  THE  RECEIVED  ACOUSTIC  PRESSURE 

A  complete  solution  to  the  problem  of  acoustic  transmission  through 
a  rough  surface  requires  that  the  probability  distribution  of  the 
received  signal  be  known  as  well  as  the  average  field  and  power.  In 
Chapter  Seven  of  his  monograph,  Beckmann  derives  the  probability  distri- 
bution of  the  sum  of  arbitrary  coplanar  random  vectors.  He  then  links 
up  this  distribution  with  the  mean  field  and  power  scattered  by  a  rough 
surface.  This  approach  provides  a  general  solution,  not  restricted  to 
mean  values,  but  giving  the  entire  distribution  about  these  mean  values. 
The  results  of  Beckmann's  analysis  can  be  applied  to  the  signal  trans- 
mitted through  a  rough  surface  and  the  same  conclusions  can  be  made 
about  both  cases. 

A  qualitative  summary  of  the  characteristics  of  this  distribution 
can  be  given  easily.  For  the  case  of  low  roughness,  R  <  1,  the  signal 
observed  in  Snell's  direction  will  have  a  near  Gaussian  distribution  of 
amplitudes.  This  result  reflects  the  addition  of  a  large  coherent 
component  to  a  small  incoherent  component.  As  the  roughness  increases 

the  amplitude  distribution  in  Snell's  direction,  K   =  0  ,  approaches 

xy 

Rayleigh,  indicating  the  absence  of  the  coherent  component.  In  the  case 
where  the  signal  is  observed  from  a  direction  away  from  Snell's  direction 
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and,  therefore,  away  from  the  predominant  smooth  surface  field,  the 
incoherent  component  dominates  and  the  distribution  will  be  near  Rayleigh 
for  any  roughness. 
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III.  THE  MODEL  EXPERIMENT 

The  theory  developed  in  the  previous  section  has  been  tested  in  a 
model  experiment  that  has  been  scaled  to  be  representative  of  the  condi- 
tions at  sea  for  the  lower  audio  frequencies.  The  significant  ratios 
a/Xp  and  L/X2  (root-mean-square  wave  height/acoustic  wavelength  and 
wave  correlation  distance/acoustic  wavelength)  take  on  values  in  the 
model  that  are  comparable  to  those  at  sea  for  the  range  of  acoustic 
frequencies  of  particular  interest.  A  result  of  this  scaling  has  been 
the  choice  of  L/A?  ratios  that  do  not  satisfy  the  requirements  of  the 
Kirchhoff  approximation.  This,  of  course,  should  be  the  objective  of  a 
model  experiment;  to  study  cases  that  may  be  beyond  the  capability  of 
theorists,  but  of  great  practical  importance. 

A.   DESIGN  OF  THE  MODEL  -  SCALING  CONSIDERATIONS 

In  order  to  design  the  particular  configuration  to  be  used  in  the 
model  experiment  some  compromises  were  made  between  what  was  considered 
ideal  and  what  was  physically  realizable  with  the  equipment  and  facilities 
available.  The  fundamental  requirement  was  to  achieve  the  desired  a/  \^ 
and  L/x  ratios.  After  that  it  was  attempted  to  satisfy  as  nearly  as 
possible  the  assumptions  made  in  deriving  the  theoretical  results.  The 
basic  configuration  of  the  experiment  consisted  of  a  source  of  airborne 
acoustic  waves  located  above  an  anechoic  tank  and  transmitting  through 
an  aperture  so  as  to  create  a  limited  ensonified  area  on  the  water  sur- 
face. A  model  sea  surface  was  generated  by  directing  a  wind  stream  over 
the  surface  of  the  anechoic  tank.  The  transmitted  acoustic  signal  was 
received  by  a  hydrophone  located  far  enough  from  the  ensonified  area  to 
be  considered  in  the  far  field. 
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The  Underwater  Acoustics  Laboratory  at  the  Naval  Postgraduate  School 
is  equipped  with  an  anechoic  tank  configured  for  generating  a  model  sea 
surface.  The  development  of  this  facility  has  taken  place  over  a  number 
of  years  and  the  model  sea  surfaces  produced  have  proven  to  be  quite 
similar  to  a  well-developed  sea  of  gravity  waves.  It  was  decided  to 
scale  the  acoustic  frequency  to  obtain  the  desired  ratios  of  a/A2 
and  L/Ap  rather  than  attempt  to  modify  the  proven  Sea  Exciter, 
(SEX  MOD.  2).  The  model  sea  surfaces  available  range  from  an  rms  wave 
height  of  about  0.1  cm  to  0.5  cm.  The  correlation  lengths  obtained  are 
scaled  to  those  of  ocean  waves  in  approximately  the  same  ratio  as  the 
rms  wave  heights,  provided  seas  with  swell  (low-frequency  gravity  waves 
propagated  from  long  distances)  are  not  considered.  In  order  to  choose 
the  range  of  acoustic  frequencies  to  be  used  the  a/Ap  ratio  was  chosen 
as  the  criterion  and  the  Pierson-Moskowitz  [4]  Ocean  Wave  Spectrum  was 
integrated  to  get  representative  rms  wave  heights  for  various  wind  speeds 
At  12  knots,  the  rms  wave  height  at  sea  should  be  about  19  cm.;  using 
this  figure  for  a  and  a  low  frequency  limit  on  the  order  of  100  Hz  (A?  = 
1500.  cm)  a  a/Ap  ratio  of  at  least  as  small  as  19/1500  should  be 
achieved  in  the  model.  Frequencies  at  least  as  high  as  100kHz  are  of 
interest  at  sea  and  for  the  same  sea  the  ratio  a/Ap  =  19/150  .  Since 
a  model  sea  with  rms  wave  height  of  .19  cm  can  be  achieved,  it  can  be 
seen  that  the  scale  factor  should  be  about  100  and  the  acoustic  frequen- 
cies used  in  the  model  should  range  from  about  10  kHz  to  100  kHz. 

B.  THE  ACOUSTIC  ARRANGEMENT 
1 .  The  Anechoic  Tank 

The  model  experiment  was  carried  out  in  an  anechoic  tank  24  feet 
long,  6  feet  wide  and  8  feet  deep.  The  tank  was  designed  to  be  anechoic 
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at  frequencies  above  20  kHz  by  use  of  a  lining  of  sound  absorbent 
aluminum-loaded  rubber.  Since  it  was  desired  to  use  frequencies  as  low 
as  10  kHz,  it  became  apparent  that  the  anechoic  qualities  of  the  tank 
would  have  to  be  studied  and  improved.  An  experiment  to  evaluate  the 
qualities  of  the  anechoic  treatment  was  performed  to  determine  how  good 
(or  bad)  it  was  at  frequencies  below  20  kHz.  This  experiment  involved 
transmitting  tone  burst  with  an  electrostatic  transducer  at  normal 
incidence  to  the  tank  bottom.  An  Atlantic  Research  LC-32  hydrophone 
located  between  the  transmitter  and  the  bottom  was  used  to  determine  the 
ratio  of  the  bottom  reflected  pulse  to  the  outgoing  pulse.  Then,  main- 
taining the  same  relative  geometry  but  inverting  the  transducers,  the 
transmitted  signal  was  directed  upward  past  the  hydrophone  and  the  ratio 
of  the  reflected  pulse  from  the  smooth,  free  surface  to  the  outgoing 
pulse  was  determined.  From  these  data  an  estimate  of  the  normal  inci- 
dence reflection  loss  for  the  tank  bottom  was  computed  and  the  results 
are  given  in  Table  2. 

A  second  experiment  was  conducted  to  study  the  reverberant  field 
in  the  tank.  A  continuous  wave,  10  kHz  tone  was  transmitted  by  a  loud- 
speaker above  the  tank  with  the  surface  acoustically  masked  except  for 
a  25  cm  square  area.  The  incident,  nearly  plane,  acoustic  wave  on  the 
smooth  water  surface  was  at  normal  incidence  and  the  field  in  the  water 
was  probed  vertically  beneath  the  ensonified  area  with  an  LC-32  hydro- 
phone. This  yielded  a  plot  that  showed  a  spherically  spreading  field 
with  a  small  but  not  insignificant  standing  wave  that  became  larger  close 
to  the  tank  bottom.  Figure  2  shows  this  plot. 

A  study  by  Lastinger  [5]  indicates  that  water-saturated  wood  can 
be  an  effective  absorber  of  low  frequency  (3-8  kHz)  sound.  In  particular 
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TABLE  2 

Tabulation  of  normal  incidence  reflection  coefficient  for  the 
aluminum-loaded  rubber  anechoic  lining  installed  on  the  tank  bottom, 


Note:  p     reflected  pressure  amplitude 


p.     incident  pressure  amplitude 
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he  reports  that  redwood  can  provide  an  attenuation  approaching  35  dB/meter 
at  8  kHz  for  sound  propagated  at  45°  to  the  grain.  For  sound  propagated 
at  0°  to  the  grain  the  acoustic  characteristics  are  said  to  be  similar 
to  those  obtained  for  45°  orientation.  Perpendicular  propagation  pro- 
vided the  least  absorption. 

Since  the  model  experiment  would  result  in  sound  directed  primarily 
at  the  tank  bottom,  it  was  felt  that  additional  treatment  of  the  bottom 
was  necessary  and  that  a  layer  of  redwood  between  the  concrete  bottom  and 
the  absorbent  rubber  lining  would  be  useful.  Redwood  boards  8  by  4  by  1 
inch  were  cut  and  wet  by  first  being  placed  under  water  in  a  large, 
sealed,  metal  container  that  was  then  evacuated  to  a  nominal  vacuum  of 
29  inches  of  Hg  for  approximately  72  hours.  These  boards  were  then  wired 
together  in  blocks  8  by  8  by  4  inches  and  placed  under  water  in  the  four 
foot  deep  runways  between  the  two  anechoic  tanks  in  the  Acoustics  Lab. 
When  enough  saturated  redwood  blocks  had  been  fabricated  to  cover  approxi- 
mately 12  feet  of  the  tank  bottom,  the  tanks  were  drained  and  the  redwood 
installed  with  the  grain  oriented  vertically.  This  provided  a  4  inch 
deep  layer  of  redwood  on  12  feet  by  6  feet  of  the  bottom  beneath  the  area 
that  would  be  ensonified  during  the  model  experiment.  Based  on  Lastinger's 
study  this  should  have  provided  additional  attenuation  of  7  dB.  The 
absorbent  rubber  lining  was  placed  on  top  of  the  redwood  blocks,  which 
reduced  the  depth  of  the  tank  to  approximately  190  cm.  After  the  redwood 
was  in  place  a  reverberation  experiment  identical  to  that  described  above 
was  conducted  and  a  noticeable  improvement  was  evident.  Figure  2  is  a 
plot  comparing  the  field  in  the  anechoic  tank  before  and  after  the  instal- 
lation of  redwood. 

The  acoustic  data  in  the  model  experiment  was  taken  approximately 
two  months  after  the  redwood  was  in  place  and  at  that  time  an  additional 
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test  was  made  to  evaluate  the  significance  of  the  backscattered  signal 
from  the  tank  bottom.  An  electrostatic  transducer  with  an  active  area 
approximately  equal  to  that  of  the  ensonified  area  used  in  the  model 
experiment  was  located  just  below  the  surface  in  the  same  position  as 
the  ensonified  area.  The  receiving  hydrophone  was  positioned  as  it  had 
been  during  the  experiment.  The  arrangement  approximated  the  conditions 
that  were  present  during  the  experiment  but  permitted  the  utilization 
of  a  source  in  a  pulsed  mode  in  order  to  observe  the  magnitude  of  the 
backscattered  signal  relative  to  outgoing  signal.  Measurements  were 
taken  for  0°,  30°  and  60°  incidence  with  the  bottom  and  the  outgoing 
pulse  amplitude  received  at  the  on-axis  hydrophone  was  compared  to  the 
backscattered  pulse  amplitudes.  Table  3  is  a  summary  of  the  results 
of  this  test.  The  results  for  10  kHz  and  the  normal  incidence  can  be 
compared  with  the  first  evaluation  of  the  loss  at  the  bottom  outlined 
above  and  an  estimate  of  the  attenuation  of  the  redwood  can  be  made. 

The  20  dB  difference  between  outgoing  and  backscattered  pulses 
can  be  attributed  to  spreading  losses,  losses  introduced  by  the  metal - 
loaded  rubber  and  losses  introduced  by  the  redwood.  Since  the  hydrophone 
was  located  about  90  cm  below  the  source  and  100  cm  above  the  concrete 
tank  bottom,  the  difference  due  to  spreading  can  be  estimated  to  be  about 
10  dB.  The  loss  introduced  by  the  rubber  anechoic  lining,  when  rigidly 
backed,  has  been  evaluated  as  5  dB.  If  these  figures  are  used,  then  it 
appears  that  the  redwood  treatment  has  resulted  in  a  5  dB  improvement  at 
10  kHz.  A  similar  analysis  for  20  kHz  indicates  that  the  addition  of 
redwood  has  apparently  reduced  the  effectiveness  of  the  anechoic  treatment 
by  7  dB.  A  possible  explanation  may  be  that  the  rubber  absorbant  material 
is  less  effective  when  not  rigidly  backed.  Although  it  is  difficult  to 
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TABLE  3 


Tabulation  of  backscattered  amplitudes  relative  to  outgoing 
amplitudes  received  by  hydrophone  at  one  meter  and  e?  degrees  from 
center  of  ensonified  area.  Redwood  layer  has  been  installed  on  tank 
bottom  with  anechoic  lining  on  the  redwood.  The  first  backscattered 
return  was  from  the  bottom  —  the  second  was  predominantly  from  the 
water  surface. 


FREQUENCY 

e2 

FIRST  RETURN 

SECOND  RETURN 

(kHz) 

10.0 

0 

0.10 

0.07 

13.5 

0  . 

0.10 

0.06 

17.5 

0 

0.09 

0.02 

20.0 

0 

0.07 

0.02 

24.0 

0 

0.06 

0.02 

30.0 

0 

0.05 

0.03 

40.0 

0 

0.04 

0.02 

10.0 

30 

0.10 

0.10 

13.5 

30 

0.06 

0.04 

17.5 

30 

0.06 

0.03 

20.0 

30 

0.04 

0.01 

24.0 

30 

0.03 

0.01 

30.0 

30 

0.01 

0.01 

40.0 

30 

0.01 

0.01 

10.0 

60 

0.09 

0.09 

13.5 

60 

0.04 

0.03 

17.5 

60 

0.02 

0.01 

20.0 

60 

0.01 

0.01 

24.0 

60 

0.01 

0.01 

30.0 

60 

0.01 

0.01 

40.0 

60 

0.01 

0.01 
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estimate  the  nature  of  the  reverberant  field  that  existed  during  the 
rough  surface  transmission  experiments,  a  conservative  evaluation  of  this 
problem  would  seem  to  indicate  that  the  data,  particularly  the  mean 
square  acoustic  pressure  data,  should  be  subject  to  close  examination.  A 
more  detailed  discussion  of  this  problem  will  be  given  in  a  later  section. 
2.  The  Ensonified  Area 

The  theory  is  expressed  in  terms  of  a  limited  area  that  has  a 
uniform  acoustic  field  over  dimensions  2X  by  2Y  and  zero  acoustic  field 
outside  these  boundaries.  In  addition  the  incident  wave  is  required  to 
be  plane  and  the  area  ensonified  must  be  larger  than  the  acoustic  wave- 
length in  water.  The  assumption  that  the  receiver  is  in  the  far  field 
necessarily  limits  the  size  of  the  ensonified  area,  while  the  description 
of  the  rough  surface  requires  a  Gaussian  distribution  of  heights  and  an 
isotropic  surface  that  has  a  Gaussian  spatial  correlation  function. 
Other  than  the  desire  to  test  the  theory  for  a  realistic  sea  surface  that 
is  anisotropic  and  does  not  have  the  prescribed  spatial  correlation,  it 
was  felt  that  the  experimental  ensonified  area  should  be  a  reasonable 
approximation  to  the  idealized  description  of  the  theory  and  yet  flexible 
enough  in  design  to  permit  the  investigation  of  some  of  the  theoretical 
limitations. 

A  number  of  ideas  for  implementing  this  critical  component  of 
the  experiment  were  tried  out  and  rejected.  A  parabolic  reflector  20 
inches  in  diameter  with  an  acoustic  feedhorn  at  the  focus  failed  to  pro- 
duce a  satisfactory  field  primarily  because  the  horn  required  to  give 
sufficient  directivity  was  so  large  as  to  severely  disturb  the  reflected 
wave.  This  approach  was  later  successful  when  a  reflector  from  a  portable 
microwave  relay  station  was  obtained  and  used  in  conjunction  with  the 
basic  arrangement.  The  result  of  this  combination  will  be  discussed  later, 
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Another  source  for  generating  a  limited  ensonified  area  was  an 
8  foot  long  fiberglass-lined  plywood  frame  enclosure  with  the  Altec  740A 
driver  and  catenoidal  horn  mounted  at  one  end  and  directed  through  the 
tube  of  fiberglass  absorbing  material  designed  to  shape  the  desired 
field.  An  arrangement  in  the  anechoic  chamber  that  tracked  a  microphone 
across  the  mouth  of  the  enclosure  permitted  the  measurement  of  the 
amplitude  and  relative  phase  of  the  resultant  field.  The  performance  of 
this  enclosure  was  improved  by  flaring  the  lower  end  so  as  to  form  a 
sort  of  fiberglass  horn  20  by  20  inches  at  one  end  and  48  by  48  inches 
at  the  other.  The  bottom  end  had  a  wire  grid  supporting  more  fiberglass 
absorbing  material  to  form  an  aperture  that  could  be  adjusted  to  any 
desired  size.  The  enclosure  was  installed  above  the  center  of  the 
anechoic  tank  supported  by  9/16  inch  diameter  steel  rods  pivoting  from 
two  steel  angle  A-frames  along  either  side  of  the  tank.  Although  this 
system  seemed  reasonably  workable,  it  was  decided  to  go  to  a  similar 
arrangement  that  provided  more  positive  masking  of  the  surface  outside 
the  desired  region  of  ensonification. 

a.  The  Basic  Arrangement 

The  basic  arrangement  finally  settled  on  consisted  of  a  56 
by  56  inch  masking  pad  placed  over  the  anechoic  tank  with  an  aperture  in 
the  center  that  was  50  cm  square  but  fitted  with  an  insert  that  could 
reduce  it  to  25  cm  square.  The  pad  was  constructed  of  three  inches  of 
fiberglass  absorbing  material  on  either  side  of  a  1/4  inch  plywood 
center.  The  underneath  side  of  the  pad  was  covered  with  a  protective 
sheet  of  1  mil  mylar.  The  fiberglass  sound  absorbing  material  was  tested 
in  the  anechoic  chamber  and  found  to  introduce  a  transmission  loss  of 
about  7  dB  per  inch  for  the  10-20  kHz  frequency  range  and  no  detectable 
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reflective  properties.     The  total  transmission  loss  of  the  pad  including 
the  plywood  septum  was  approximately  65  dB. 

The  masking  pad  was  installed  6  inches  above  the  water  level 
in  the  tank  so  as  to  permit  the  sea  generating  wind  stream  to  pass  under 
it.     Ten  feet  above  the  pad  was  the  Altec  driver  and  horn  mounted  on  a 
three  foot  long  fiberglass-lined  enclosure  that  was  pivoted  from  the 
A-frame  supports  described  previously.     The  height  of  the  mouth  of  the 
catenoidal   horn  was  approximately  305  cm  above  the  water  level.     A  light 
plywood  frame  strung  with  a  three  inch  square  grid  of  fine  wire  that 
supported  fiberglass  absorbing  material  was  positioned  over  the  tank 
just  beneath  the  source  enclosure  and  formed  the  roof  of  a  fiberglass 
enclosure  that  was  fitted  over  the  masking  pad. 

The  A-frames  and  therefore  the  acoustic  source  could  be 
positioned  along  the  length  of  the  tank  to  allow  various  angles  of  inci- 
dence from  the  source  to  the  ensonified  area.     The  aperture  insert  could 
be  oriented  perpendicular  to  the  indident  field  for  angles  within  about 
20  degrees  of  the  vertical. 

An  LC-32  hydrophone,  mounted  on  a  traverse,  was  swept  across 
the  ensonified  area  to  measure  a  cross-section  of  the  amplitude  and 
relative  phase  of  the  incident  field.     The  traverse  could  be  positioned 
to  sweep  the  hydrophone  just  below  the  water  surface   (~1   cm)   and  either 
along  the  longitudinal   axis  or  across  the  anechoic  tank.     A  potentiometer 
linked  to  the  drive  of  the  traverse  was  used  to  supply  position  informa- 
tion to  the  x-axis  of  a  Varian  F-100  X-Y  Plotter.     The  y-axis  was  supplied 
with  a  signal   from  either  a  Hewlett  Packard  (Mod.   7035)  Logarithmic  Con- 
verter or  an  Acton  Laboratories  Phase  Meter  (Mod.  AB)   in  order  to  plot 
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either  the  amplitude  or  relative  phase  distribution  across  the  ensoni- 
fied  area.  Figure  3  is  a  sketch  of  the  basic  acoustic  arrangement. 

The  properties  of  the  field  produced  on  the  surface  were 
investigated  with  the  arrangement  described  above.  In  general,  it  was 
found  that  the  ensonified  area  for  all  the  frequencies  used  was  charac- 
terized by  a  central  area  of  nearly  uniform  sound  level  that  abruptly 
decreased  with  a  slope  of  about  2  1/2  dB  per  centimeter  as  the  probe 
was  moved  beyond  the  limits  of  the  aperture.  The  dimensions  of  the 
ensonified  area,  taken  at  the  3  dB  points  was  approximately  24  cm  by 
24  cm.  The  Acton  Laboratories  Phase  Meter  compared  the  phase  of  the 
driving  voltage  to  the  source  with  that  of  the  hydrophone  and  provided 
a  d.c.  voltage  proportional  to  the  relative  phase  that  was  plotted  on 
the  x-y  recorder  as  the  hydrophone  was  traversed  across  the  ensonified 
area.  The  results  showed  wavefront  curvatures  that  were  consistent  with 
a  spherically  diverging  wave  from  a  source  at  a  distance  on  the  order  of 
three  meters. 

A  calculation  of  the  expected  wavefront  curvatures  can  be 
made  using  the  formula  for  the  sagitta  of  an  arc: 

(2R  -  x)  x  =  d2 

where  R  is  the  radius  of  curvature,  2d  is  chord  of  the  arc  and  x 
is  the  sagitta  of  the  arc.  Using  R  =  300  and  d  =  12  ,  the  sagitta 
turns  out  to  be  .24  cm.  If  the  curvature  of  the  incident  wavefront  is 
expressed  in  terms  of  the  fraction  of  a  wavelength  phase  difference 
between  the  center  of  the  ensonified  area  and  the  edge,  assuming  normal 
incidence,  then  the  ratio  of  the  sagitta  of  the  arc  and  the  acoustic 
wavelength  should  give  this  fraction.  For  an  incident  10  kHz  wave, 
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therefore,  the  expected  value  of  wavefront  curvature  is  about  .07  A,  . 
This  value  compares  to  the  0.1  X,  value  obtained  in  the  experiment  and 
plotted  in  Figure  4.  The  highest  frequency  used  with  this  source  and 
configuration  was  24.3  kHz.  The  wavefront  curvature  obtained  at  this 
frequency  was  approximately  .17  A,  .  A  criterion  for  permitting  the 
plane  wave  approximation  that  has  been  given  by  Medwin  and  Clay  [6]  is 
a-maximum  curvature  of  A/8  .  In  this  experiment  this  has  been  exceeded 
somewhat  and  an  opportunity  has  been  provided  to  test  the  criterion. 

An  indicator  of  how  closely  the  experimental  ensonified  area 
approximated  the  theoretical  description  is  given  by  a  comparison  of  the 
smooth  surface  radiated  field  measured  in  the  experiment  with  that  pre- 
dicted by  the  theory.  Since  for  a  smooth  surface  the  theoretical 
predictions  reduce  to  a  relatively  straight-forward  expression  for  the 
radiation  pattern  of  a  piston  with  a  particular  linear  phase  shading, 
there  should  be  reasonably  good  agreement  at  least  for  the  main  features. 
The  plots  in  Figures  5-8  are  representative  of  those  obtained  with  the 
basic  acoustic  arrangement.  Each  plot  compares  the  theoretical  and 
experimental  directivity  patterns  of  the  ensonified  area  for  a  given 
frequency  (10.0  kHz  or  24.3  kHz)  and  angle  of  incidence,  9,,  (0°,  6°, 
12°  or  16°).  The  angles  are  those  used  in  the  rough  surface  study  and 
were  selected  to  be  below  and  above  the  critical  angle,  9  ,  for  infinite 
plane  waves  incident  on  a  plane  surface  (about  13°).  The  incident 
angles  were  measured  by  the  geometry  of  the  arrangement  and  are  estimated 
to  be  accurate  to  within  one  degree.  Since  angular  position  of  the  main 
lobe  is  very  sensitive  to  small  changes  of  incident  angle,  some  of  the 
theoretical  patterns  have  been  rotated  as  much  as  5°  to  obtain  the  best 
fit. 
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FIGURE  4.  PHASE  AND  AMPLITUDE  CROSS-SECTION  AT  ENSONIFIED  AREA 
SOURCE:  Altec  driver  and  horn  at  10.0  kHz 
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b.  The  Parabolic  Reflector 

In  an  effort  to  extend  the  frequency  range  of  the  experiment 
it  became  necessary  to  use  a  new  acoustic  source.  Two  requirements  were 
placed  on  the  new  source;  it  had  to  be  able  to  transmit  sufficient  sound 
pressure  levels  at  the  higher  frequencies  and  it  should  have  produced  a 
nearly  plane  wave  incident  on  the  ensonified  area.  A  parabolic  section 
reflector  from  a  AN/TRC-27  Portable  Microwave  Relay  Set  was  obtained 
which,  in  conjunction  with  a  Western  Electric  640AA  condenser  micro- 
phone used  as  a  transmitter,  produced  a  most  satisfactory  incident  wave 
and  permitted  the  extension  of  the  experiment  to  43  kHz.  The  microwave 
horn  was  removed  from  its  mount  at  the  focal  point  of  the  reflector  and 
the  640AA  transmitter  was  installed.  The  reflector  was  positioned 
approximately  1.8  meters  above  the  water  surface  and  directly  over  the 
aperture  in  the  masking  pad.  Provisions  for  making  fine  adjustments  of 
the  reflector  unit  were  built  into  the  reflector.  Adjustment  of  the 
reflector  until  a  maximum  signal  was  obtained  on  the  receiving  hydro- 
phone located  at  62  =  0°  was  easily  accomplished.  The  surface 
hydrophone  was  traversed  across  the  ensonified  area  and  an  amplitude  and 
phase  cross-section  was  obtained.  The  amplitude  distributions  were 
essentially  the  same  as  had  been  obtained  with  the  original  source  but 
the  phase  distributions  showed  no  discernable  curvature  even  at  37  kHz. 
See  Figure  9  for  the  plots  obtained  at  37  kHz. 

The  source  level  of  the  W.E.  640AA  transmitter  was  on  the 
order  of  0  dB  re  1  ybar  at  30  cm  when  polarized  at  180  volts  and  driven 
with  30  vrms.  The  signal  levels  obtained  at  the  hydrophone  were  consid- 
erably lower  than  those  experienced  with  the  Altec  740A  source.  Since 
the  signal  levels  were  about  20-25  dB  below  those  obtained  previously, 
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FIGURE  9.  PHASE  AND  AMPLITUDE  CROSS-SECTION  OF 
"PLANE  WAVE"  AT  ENSONIFIED  AREA 

SOURCE:  Parabolic  Reflector  and 
W.E.  640AA  at  37.0  kHz 
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the  noise  in  the  receiving  electronics  precluded  computation  of  the 
relative  mean  square  pressure  until   the  narrower  filter  of  a  wave 
analyzer  was  inserted  to  improve  the  signal-to-noise  ratio.     This  un- 
corrected additive  noise  was  not  a  problem  in  computing  the  relative 
mean  acoustic  pressure  but  the  low  signal   levels  did  cause  difficulties 
because  of  a  small,  pervasive,  correlated  electrical   feedthrough  from 
the  transmitting  electronics  that  became  a  factor  when  the  acoustic 
level   decreased. 

c.     The  Point  Source 

In  a  supplementary  experiment  the  transmission  of  a  nearly 
omnidirectional   point  source  through  the  rough  air-sea  interface  without 
an  unlimited  ensonified  area  was  studied.     The  source  used  was  the 
Western  Electric  640AA  condenser  transmitter,  which  has  an  active  element 
diameter  of  1.7  cm.     Reference  [7]  discusses  the  characteristics  of  this 
microphone  used  as  a  transmitter.     The  field  of  the  transmitter  is 
piston-like  and  for  the  frequencies  used  would  be  a  reasonable  approxi- 
mation to  an  omnidirectional   point  source  out  to  angles  greater  than  the 
critical   angle  defined  by  Snell's  Law.     Figures  10  and  11   show  the  free 
field  patterns  of  the  transmitter  used  obtained  at  10  and  20  kHz  in  the 
anechoic  chamber. 

In  this  experiment  the  masking  pad  was  removed  and  the  micro- 
phone was  positioned  1  meter  above  the  water  surface  and  over  the  center 
of  the  tank.     An  enclosure  of  fiberglass  absorbing  material  was  constructed 
over  the  tank  with  the  sides  extending  down  to  within  about  3  cm  of  the 
water  surface  and  the  top  supported  by  a  fram  with  a  fine  wire  grid  just 
above  the  height  of  the  transmitting  microphone.     This  left  about  16 
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FIGURE  10. 


FREE  FIELD  DIRECTIVITY  PATTERN  AT  1  METER 
FROM  W.E.  640AA  TRANSMITTER  AT  10.0  kHz 
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FIGURE  11. 


FREE  FIELD  DIRECTIVITY  PATTERN  AT  1  METER 
FROM  W.E.  640AA  TRANSMITTER  AT  20.0  kHz 
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square  feet  of  water  surface  exposed  to  the  nearly  omnidirectional 
source  and  reduced  the  contributions  of  secondary  scattering  to  the 
field  over  the  ensonified  area. 

3.  The  Receiver  and  Associated  Electronics 

An  Atlantic  Research  LC-32  hydrophone  with  active  element  dimen- 
sions of  3/4  by  1/4  inches  was  used  as  the  reciever.  The  hydrophone  was 
mounted  on  a  L-shaped  support  that  pivoted  at  the  water  surface  and 
positioned  the  hydrophone  in  the  plane  passing  through  the  source,  the 
center  of  the  ensonified  area,  and  the  central  axis  of  the  tank.  A 
potentiometer  pickoff  at  the  pivot  point  supplied  a  d.c.  voltage  propor- 
tional to  the  angle  between  the  center  of  the  ensonified  area  and  the 
hydrophone.  This  arrangement  was  used  to  swing  the  hydrophone  through 
a  180°  arc  to  obtain  the  smooth  surface  polar  plots  and  also  to  position 
the  hydrophone  at  a  given  angle  for  measuring  the  relative  average 
transmitted  acoustic  field  and  intensity  for  the  rough  surface. 

An  NUS  low  noise  preamplifier  (Model  #2110)  amplified  the  hydro- 
phone signal  by  30  dB.  Since  the  sound  pressure  levels  in  the  experiment 
were  low  (on  the  order  of  20  dB  re  1  ybar)  and  the  hydrophone  sensitivity 
was  about  -110  dB  re  1  volt/ybar,  an  additional  60  dB  of  amplification  by 
a  Hewlett-Packard  465A  Precision  Amplifier  was  required  to  bring  the 
signal  level  to  the  one  volt  range.  An  SKL  variable  electronic  filter 
was- used  to  bandpass  filter  the  signal  with  an  effective  bandwidth  of 
about  1  kHz.  A  variation  on  the  above  filtering  technique  was  the  use  of 
a  Hewlett-Packard  3590A  wave  analyzer  with  a  bandwidth  of  100  Hz.  This 
was  used  to  obtain  the  mean  square  accoustic  pressures  when  the  Western 
Electric  640AA  transmitter  was  used.  The  filtered  output  was  supplied 
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to  one  input  of  the  PAR  Model  101  Correlation  Function  Computer  which 
was  used  to  process  all  of  the  acoustic  data. 

The  requirement  that  the  receiver  be  in  the  far  field  was  esti- 
mated by  applying  the  far  field  criterion  for  a  piston  source;  i.e., 

p 

the  far  field  was  taken  to  be  at  a  range  greater  than  2a  A2  from  the 

source  with  half  width  "a".  For  the  highest  frequency  used,  43  kHz, 
2a  /A2  is  approximately  84  cm.  Accordingly  the  hydrophone  range  was 
chosen  to  be  at  1  meter;  a  position  in  the  far  field  and  also  not  too 
close  to  the  more  reverberant  field  near  the  bottom  of  the  anechoic  tank. 

C.  THE  MODEL  SEA  SURFACE 

The  north  anechoic  tank  in  the  Underwater  Acoustics  Laboratory  is 
equipped  to  generate  a  wind-driven  model  sea  surface.  The  wind  generator 
is  located  at  the  east  end  of  the  tank  and  consists  of  five  separately 
controllable  centrifugal  fans  feeding  a  common  plenum  that  discharges  a 
uniform  wind  stream  across  the  width  of  the  tank.  The  output  of  the  wind 
generator  is  confined  to  the  water  surface  by  plexiglass  covers  suspended 
six  inches  above  the  water  level.  The  sea  surface  is  effectively  termi- 
nated at  the  west  end  of  the  tank  by  a  sloping  beach  made  of  redwood 
chips  packed  loosely  in  a  plastic  insect  screen 

The  author  assisted  Dr.  W.  M.  Wright  in  a  study  of  two  different  sea 
surfaces  that  were  later  to  be  used  in  this  model  experiment.  The  wave 
probe  system  used  was  of  a  design  similar  to  that  used  by  Mayo  [8]  and 
had  been  fabricated  and  tested  by  Dr.  Wright.  The  two  surfaces,  one  a 
"high  sea"  and  the  other  a  "low  sea",  were  studied  in  some  detail  and 
spatial  correlation  functions,  directional  surface  wave  spectra,  and  wave 
height  probabilities,  as  well  as  the  mean  square  wave  heights,  were 
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obtained  for  a  location  centered  at  the  ensonified  area.  Later  in  the 
experiment,  when  it  became  desirable  to  generate  other  sea  surfaces  with 
particular  characteristics,  the  author  used  the  same  wave  probe  system 
to  make  a  limited  study  of  these  supplementary  seas. 
1 .  The  Wave  Probe  System 

Reference  [9]  gives  a  detailed  discussion  of  the  design,  con- 
struction, and  testing  of  the  wave  probes  and  associated  electronics. 
The  probes  were  based  on  the  expectation  that  the  conductance  between 
two,  long  conducting  cylinders  is  proportional  to  the  depth  of  immersion 
of  the  cylinders.  The  circuit  incorporating  the  wave  probe  was  designed 
to  detect  the  water  level  by  measurement  of  the  conductance  between  the 
probe  electrodes.  Figure  11  shows  the  complete  circuit  for  wave  probe 
"A".  The  probe  is  driven  by  an  oscillator  with  low  source  resistance 
and  the  current  is  measured  by  means  of  the  10  ohm  resistor.  A  passive 
filter  removes  any  contribution  at  the  carrier  frequency  of  wave  probe 
"B"  from  the  amplitude-modulated  signal.  This  signal  is  amplified, 
detected  and  filtered  to  produce  the  wave  record  which  is  again  amplified 
before  processing  by  the  PAR  Correlator.  The  "AC  Input"  of  the  corre- 
lator was  used,  which  provided  a  lower  cutoff  frequency  of  about  0.2  Hz. 
Data  recorded  on  the  PI-6200  FM  tape  recorder  was  passed  through  a  50  yF 
blocking  capacitor  into  the  50  kft  input  impedance  of  the  tape  recorder. 

When  it  was  desired  to  make  calibrated  measurements  with  the  wave 
probe,  a  procedure  was  followed  that  provided  a  calibration  of  the  entire 
wave  probe  system  including  the  correlator.  This  was  necessary  since  it 
was  apparent  that  the  sensitivity  of  the  probe  was  not  stable  over  a  long 
term;  possibly  due  to  a  changing  conductivity  in  the  water  in  the  anechoic 
tank. 
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The  calibration  procedure  consisted  of  oscillating  the  probe 
vertically  in  the  anechoic  tank  with  a  motor  driven  apparatus.  A  sinu- 
soidal variation  of  immersion  with  time  was  obtained,  with  a  frequency 
of  2.72  Hz  and  a  peak-to-peak  amplitude  of  2  cm.  The  signal  generated 
by  this  method  was  applied  to  the  correlator  and  the  autocorrelation  at 
x  =  0,  the  mean  square  voltage,  was  computed.  The  resulting  system 
calibration  was  given  by  .707  cm/[C(0)]  '       volt  and  averaged  about 
.55  cm/volt.  The  d.c.  level  of  the  wave  record  was  noted  and  the 
correct  immersion  depth  was  obtained  after  calibration  by  adjusting  the 
depth  to  produce  the  same  d.c.  voltage.  Calibration  for  the  data  to  be 
processed  by  the  digital  computer  was  accomplished  by  recording  a  ten 
minute  record  of  the  calibration  signal  at  the  beginning  of  the  tape 
containing  the  wave  records.  Again,  the  same  d.c.  level  observed  during 
the  calibration  run  was  maintained  during  the  recording  of  the  wave 
record  by  adjusting  the  wave  probe  depth. 

2.  Statistical  Considerations  in  Processing  the  Wave  Record 

Since  the  rms  wave  height  is  a  fundamental  quantity  in  describing 
the  "roughness"  of  the  air-sea  interface  in  terms  of  the  theory,  the  meas- 
urement of  this  quantity  was  considered  fundamental  to  the  study  of  the 
sea  surface.  Studies  in  the  past  had  indicated  that  the  mean  square  wave 
height  of  the  model  seas  exhibited  temporal  fluctuations  that  appeared 
to  be  statistical  in  nature  [8]. 

Bendat  and  Piersol  [10]  give  a  discussion  of  the  statistical 
errors  to  be  encountered  in  parameter  estimates.  For  the  variance  of  the 
estimate  for  the  mean  square  of  a  zero  mean  process  they  find: 

Var  [ax2]  -  |   /  Cx2(t)  dx 
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A.      n 

for  large  T,  where  |t|  «  T  .  It  can  be  seen  that  a    is  a  consistent 

2  ~  2 

estimate  of  ax   since  Var[ax  ]  will  approach  zero  as  T  approaches 

o 
infinity  assuming  C  (t)  is  absolutely  integrable  over  (-°°,<*0  . 

If  the  wave  record  is  for  simplicity  taken  to  be  band  limited 

white  noise  with  a  bandwidth,  B,  and  length,  T  .  Then  it  can  be  shown 

that  the  above  becomes 

-  2      Rx2<°)      ax4 
Var[a  2]  -  -^ ^  JL. 

*  BT  BT 


or 


2  Var  K2J     ^     1 

=     —<  -     BT 

X 

The  square  root  of  this  last  expression  gives  the  normalized  standard 
error,  e   . 

For  the  case  of  an  RC  integrator  as  used  in  the  PAR  Correlator 
the  averaging  time  T  used  above  can  be  given  by  the  formula: 

-Tr/RC     2 


2  RC  [1  -  e     '       ] 
[1  -  e"2T^RC] 


where  T   is  the  length  of  the  record  available.  Since  there  was  no 
reason  the  record  could  not  be  as  long  as  desired,  the  effective  averaging 
time  in  the  correlator  could  be  considered  to  be  2RC  or  40  seconds.  If 
the  effective  bandwidth  of  the  wave  record  is  taken  to  be  about  3  Hz, 
then  the  normalized  standard  error  can  be  computed  to  be,  e  =  .09  , 
which  would  account  for  the  ±20%  fluctuations  in  the  mean  square  wave 
height  reported  by  Mayo  and,  in  fact,  observed  in  this  experiment  as  well. 
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In  order  to  reduce  the  error  in  the  estimate  of  the  mean  square 
wave  height  it  was  necessary  to  increase  the  averaging  time.     This  was 
accomplished  by  plotting  a  500  second  record  of  the  mean  square  wave 
height  as  computed  by  the  PAR  Correlator  and  visually  averaging  this 
record.     It  was  found  that  this  technique  consistently  produced  repeat- 
able  results.     A  quantitative  check  on  the  stationarity  of  the  random 
sea  surface  generated  was  performed  by  processing  a  wave  record  in  a 
digital   computer.     In  this  analysis  a  20  minute  wave  record  for  each  sea 
studied  was  recorded  on  the  PI-6200  FM  tape  recorder  in  analog  form. 
This  sample  was  converted  to  digital    information  by  the  A  to  D  conversion 
capability  of  the  CI5000-XDS9300  hybrid  computer  facility  and  processed 
digitally  to  obtain  mean  square  values  for  two  ten-minute  blocks  and  the 
entire  20-minute  block.     The  values  obtained  for  the  ten-minute  blocks 
were  within  1%  of  that  obtained  for  the  twenty-minute  block  and  the 
computed  mean  square  wave  heights  agreed  with  those  obtained  through  the 
use  of  the  correlator  in  the  manner  outlined  above. 
3.     Description  of  the  Sea  Surfaces  Generated 

For  the  basic  experiment  the  two  sea  surfaces  studied  were 
characterized  as  a  "low"  sea  and  a  "high"  sea.     The  root-mean-square 
wave  heights  at  the  center  of  the  tank,  where  the  ensonified  area  was  to 
be  located,  were  determined  to  be  0.17  cm  and  0.41   cm  for  the  two  seas. 
Using  the  scale  factor  of  100  and  assuming  a  Gaussian  wave  height  distri- 
bution, ocean  waves  with  a  peak-to-peak  height  (±3a)  of  about  3.5  feet 
and  8  feet  such  as  would  be  expected  in  a  fully  developed  sea  generated 
by  12  and  18  knot  winds,  were  simulated  by  the  two  model   seas.     These 
seas  were  designated  A  and  B  respectively  and  were  used  as  the  rough 
surfaces  for  most  of  the  acoustic  measurements.     Three  supplementary 
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seas  were  generated  during  the  course  of  the  acoustic  data  collection 
in  an  attempt  to  study  an  apparent  effect  introduced  by  varying  the 
correlation  length  of  the  sea  surface.  These  seas  were  designated  C,  D, 
and  E.  A  summary  of  the  basic  parameters  of  the  five  seas  used  is  found 
in  Table  4. 

a.  Mean  Square  Wave  Height 

The  mean  square  wave  height  was  obtained  for  all  seas  used. 
As  noted  above,  a  particular  technique  was  required  to  obtain  satisfactory 
values  for  the  mean  square  wave  height  using  the  PAR  Correlator  with  its 
limited  20  second  time  constant.  The  study  of  seas  A  and  B  was  accom- 
plished approximately  two  months  prior  to  the  time  the  bulk  of  the 
acoustic  data  were  taken;  however,  mean  square  wave  heights  for  these 
seas  were  measured  again  with  the  acoustic  arrangement  in  place  over  the 
anechoic  tank  and  found  to  be  unchanged. 

An  investigation  of  the  fetch-limited  variation  of  the  wind- 
driven  seas  was  made  by  measuring  the  mean  square  wave  heights  at  positions 
along  the  center  line  of  the  tank  but  30  cm  upwind  and  30  cm  downwind  of 
tank  center.  For  both  seas  the  mean  square  wave  height  was  smaller  by 
about  5%   in  the  upwind  direction  and  larger  by  the  same  amount  in  the 
downwind  direction.  For  the  25  cm  dimensions  of  the  ensonified  area 
using  the  mean  square  value  obtained  at  the  center  of  the  area  as  an 
average  value  for  the  entire  area  seems  to  be  a  valid  approximation. 

Assignment  of  an  interval  of  high  confidence  for  the  measured 
mean  square  wave  heights  is  possible  if  it  assumed  there  was  no  consistent 
error  introduced  by  the  calibration  technique.  The  normalized  standard 
error,  e  ,  for  a  single  500  second  average  of  the  correlator-computed 
mean  square  value  is  about  .025.  Using  a  conservative  criterion  that 
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TABLE 

4 

SEA  SURFACE 

CHARACTERISTICS 

a(cm) 

Lu  <cm) 

Lc  (cm) 

0.173 

2.1 

5.2 

0.410 

3.5 

11.0 

0.170 

3.2 

- 

0.320 

4.2 

- 

0.050 

0.8 

_ 

SEA  o(cm)  Lu   (cm)             LQ  (cm)  How  Generated 

A  0.173                  2.1                   5.2  Three  back  fans 

B  0.410                  3.5  11.0  All   five  fans 

C  0.170                  3.2                  -  Four  fans  -  with  deflector 

D  0.320                  4.2                  -  Five  fans  -  with  deflector 

E  0.050                  0.8                  -  One  fan 

Definitions: 

o    is  the  root-mean-square  wave  height 

L.,    is  the  smallest  separation  for  a  correlation  of  e" 
(upwind-downwind) 


u 


L    is  the  larqest  separation  for  a  correlation  of  e 
c    (crosswind) 
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assumes  the  correct  value  must  lie  within  ±3e  of  the  estimate,  a 
figure  of  ±7.5%  is  arrived  at  for  the  accuracy  of  the  mean  square  wave 
heights.  During  the  study  of  the  two  basic  sea  surfaces  many  calibra- 
tions and  subsequent  mean  square  computations  were  performed  using  both 
the  correlator  and  the  digital  computer.  Essentially  all  those  runs 
were  well  within  7.5%  of  the  mean.  It  is  felt  that  the  point  of  the 
above  exercise  is  simply  to  demonstrate  the  outer  limits  of  the  error 
that  might  be  introduced  by  these  wave  height  measurements. 

b.  Wave  Height  Probability  Distributions 

The  theory  assumes  that  the  probability  distribution  of  wave 
heights  of  a  randomly  rough  sea  will  be  Gaussian.  This  assumption  is  a 
reasonable  one  since  the  wave  height  distributions  that  have  been  ob- 
served at  sea  are  nearly  Gaussian.  Because  ocean  waves  are  characterized 
by  sharp  peaks  and  long  smooth  troughs,  the  distributions  observed  are 
usually  skewed  with  larger  positive  values  possible  but  low  negative 
values  more  likely.  This  departure  from  the  ideal  Gaussian  distribution 
is  usually  characterized  by  higher  order  moments  defining  the  "Skewness" 
and  "Kurtosis".  Wave  height  distributions  were  obtained  for  seas  A  and 
B  and  a  detailed  discussion  of  this  phase  of  the  sea  surface  study  is 
given  in  Appendix  A  along  with  plots  of  the  resultant  distributions. 
The  distributions  are  nearly  Gaussian  and  show  the  characteristic 
skewed  shape. 

c.  Spatial  Correlation  and  the  Directional  Spectrum 

A  more  complete  description  of  a  sea  surface  than  simply 
stating  the  mean  square  wave  height  is  given  by  specifying  also  the  two- 
dimensional  spatial  correlation  function  or  its  Fourier  Transform,  the 
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normalized  directional  wave  spectrum.  The  details  concerning  the  deter- 
mination of  this  description  are  given  in  Appendix  A.  In  the  theory  it 
is  assumed  that  the  sea  is  isotropic  and  has  a  Gaussian  correlation 
function.  Figures  A3-A10  in  Appendix  A  show  that  the  model  seas  were 
close  to  the  real  ocean,  that  is,  decidedly  anisotropic. 

The  parameter  that  is  defined  in  the  theory  as  the  correlation 
length,  L,  is  not  obviously  defined  for  this  type  of  surface.  Since  an 
important  theoretical  criterion  concerning  the  shape  of  the  surface  deals 
with  the  requirement  of  the  Kirchhoff  approximation  for  the  minimum 
radius  of  curvature  to  be  much  greater  than  the  acoustic  wavelength,  a 
characteristic  length  is  taken  to  be  that  distance  at  which  the  two- 
dimensional  spatial  correlation  first  falls  to  e"  .  This  length  has  been 
designated  the  experimental  correlation  length,  L  » and  it  turns  out  that 
it  is  always  located  along  the  wind  line. 

A  correlation  length  is  presented  in  Table  4  for  all  seas. 
However,  for  seas  C,  D,  E  a  complete  spatial  correlation  was  not  taken 
and  only  enough  of  the  upwind-downwind  spatial  correlation  was  measured 
to  determine  the  separation  required  for  the  normalized  cross-correlation 
to  reach  0.368.  Although  the  primary  reason  for  introducing  these  three 
supplementary  seas  was  to  observe  the  effect  that  modifying  the  surface 
"shape"  had  on  sound  transmission,  a  time  limitation  precluded  a  more 
complete  study  of  the  characteristics  of  these  surfaces. 

Seas  C  and  D  were  both  produced  by  deflecting  the  wind  up  off 
the  water  surface  at  a  distance  of  approximately  60  inches  upwind  of  the 
ensonified  area.  This  eliminated  much  of  the  high  frequency  energy  of 
the  surface  wave.  The  resulting  seas  were  somewhat  similar  to  ocean  swell 
that  is  generated  in  one  area  and  propagates  into  another.  Sea  surface  C 
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was  adjusted  to  have  the  same  mean  square  wave  height  as  A,  but  greater 
correlation  length,  by  adjusting  the  location  of  the  deflector,  by 
varying  the  height  of  the  deflector  so  that  some  wind  could  continue 
along  the  surface,  and  also  by  selecting  an  appropriate  combination  of 
fans  to  be  operated.  This  technique  produced  a  50%  longer  correlation 
length  than  was  found  with  sea  A  and  although  the  two  surfaces  had  the 
same  mean  square  wave  height,  sea  C  appeared  considerably  smoother. 
Sea  surface  D  was  the  result  of  an  attempt  to  generate  the  largest  possible 
correlation  length  and  at  the  same  time  a  moderate  mean  square  wave  height 
somewhere  between  those  of  seas  A  and  B.  This  was  accomplished  by  re- 
locating the  ensonified  area  60  inches  downwind  of  the  tank  center  to 
allow  increased  fetch  for  sea  build-up  and  by  adjustment  of  the  wind 
deflector. 

Sea  E  was  generated  to  produce  a  small  correlation  length. 
Only  the  center  fan  was  operated  and  the  resultant  sea  observed  at  tank 
center  was  small  but  not  so  small  as  to  be  in  the  capillary  wave  regime 
(defined  to  be  where  the  wave  length  is  less  than  1.73  cm). 

D.   PROCESSING  THE  RECEIVED  ACOUSTIC  PRESSURE 

All  processing  of  acoustic  data  was  done  with  a  PAR  Model  101 
Correlation  Function  Computer.  This  instrument  is  a  special  purpose 
hybrid  computer  that  computes  100  points  of  the  correlation  function  of 
two  input  voltages.  The  correlator  converts  one  input  channel  to  a 
sampled  digital  signal  and  introduces  100  discrete  delays.  The  undelayed 
input  channel  is  multiplied  with  the  output  of  each  of  the  100  delay 
steps  and  each  product  is  averaged  in  an  RC  integrator  with  a  20  second 
time  constant.  The  result  is  an  estimate  of  the  cross-correlation  of  the 


62 


two  input  signals  computed  for  100  values  of  x  ,  the  delay  time.  The 

output  of  the  correlator  was  displayed  on  an  oscilloscope  when  all  100 

points  were  to  be  observed.  If  the  autocorrelation  at  t  =  0  was 

desired,  the  output  for  this  single  point  was  displayed  on  an  X-Y 

recorder  with  a  time  base  input  to  the  x-axis. 

The  acoustic  data  are  presented  in  terms  of  relative  quantities 

referenced  to  the  smooth  surface  value.  This  procedure  is  convenient 

and  it  yields  a  result  that  can  be  directly  compared  to  theory.  A 

simple  modification  of  Equations  (15)  and  (29)  is  all  that  is  necessary 

to  express  those  results  in  terms  of  relative  quantities. 

1 .  The  Relative  Mean  Acoustic  Pressure  Squared 
(Relative  Coherent  Intensity) 

Since  [p20  F  sine  K  x  sine  Kv]  is  the  received  acoustic  pressure 
for  the  smooth  surface  for  a  given  geometry,  the  ratio  of  the  mean  re- 
ceived acoustic  pressure  for  the  roughened  surface  to  the  smooth  surface 
pressure  can  be  written: 

<p>     =        <P2>  -     e-R/2  (32) 

P20FD 

where  D  is  the  directivity  of  the  ensonified  area,  sine  K  x  sine  K  v  , 
and  <^p^>  is  defined  as  the  relative  mean  acoustic  pressure  at  the 
receiver.  The  quantity  that  is  of  more  direct  interest  is  the  coherent 
component  of  the  acoustic  intensity  which  is  simply  the  square  of  (32) 
and  may  be  written: 

l<P2>  I2      _     .-R 
"F2D 


R_-=.    kV    c,. 


<p><p>*     .       I^Z'I         =     e"R  (33) 

2     rl     ^2 
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2.  The  Relative  Mean  Square  Acoustic  Pressure 
(Relative  Total  Intensity) 

A  similar  definition  can  be  made  for  the  total  relative  acoustic 

intensity. 

/  \      R    .L2   -R       Rn    -Ky»2L2/4n 

<>)  =  e  R  +  2L   e  R    z        ^_    e  xy         (34) 

x  x  AD^       n=l   nn- 

These  definitions  put  the  theoretical  predictions  in  a  form  easily 
compared  to  experiment  and  also  limit  considerations  to  the  effects  of 
the  rough  surface  only. 

3.  Data  Collection  Procedure 

The  correlator  was  used  in  a  different  mode  for  each  of  the  two 
quantities  measured.  The  procedures  used  were  simple  but  required  a 
good  deal  of  time  to  perform.  The  time  requirement  was  a  result  of 
taking  sample  averages  of  the  correlator  output  to  obtain  each  data 
point.  The  variations  of  correlated  acoustic  signal  were  generally  less 
than  those  encountered  in  the  measurement  of  the  mean  square  wave  height 
because  the  received  sound  scatters  from  many  points  of  the  surface. 
Nevertheless  it  was  apparent  that  the  20  second  integration  time  of  the 
correlator  was  too  short  to  obtain  a  good  average  for  the  acoustic  signal 
without  additional  averaging  as  described  below. 

a.  Obtaining  the  Relative  Mean  Acoustic  Pressure  Squared 

A  cross-correlation  technique  was  used  to  isolate  the  coher- 
ent component  of  the  acoustic  signal  from  the  incoherent  component.  Since 
the  driving  voltage  for  the  source  could  be  used  to  supply  a  reference 
signal  proportional  to  the  transmitted  acoustic  pressure  and  the  amplified 
and  band-pass  filtered  signal  from  the  hydrophone  was  proportional  to  the 
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received  acoustic  pressure,  a  reference  value  for  the  smooth  surface 
could  be  obtained  by  cross-correlating  those  two  signals.  The  result 
of  this  cross-correlation  was  proportional  to  the  smooth  surface  received 
acoustic  pressure.  The  computation  performed  by  the  correlator  yielded 
100  sample  points  of  the  following  integral. 


-,     20  sec 
C0(t)  =  ^    /   p2Q  F  Dcos  aj(t-x)  cos  co(t-t')  dt 


where  C  (t)  is  an  unbiased  estimate  of  the  cross-correlation  function 
o 

and  the  integral  including  the  uncorrelated  additive  noise  has  been  set 
equal  to  zero.  Smooth  surface  cross-correlation  yielded: 

Cq(t)  =  -^ cos  [WT  -  *] 

This  function  was  displayed  on  the  oscilloscope  and  the  peak-to-peak 
value  was  adjusted  to  be  full  scale  (10  cm)  on  the  scope.  Then  without 
changing  the  driving  voltage  to  the  source  the  surface  was  roughened  and 
the  resultant  correlation  function  was  observed.  It  was  noted  that  the 
phase  of  the  correlation  function  did  not  change  but  that  the  amplitude 
was  reduced  when  the  surface  was  roughened.  In  most  cases  the  gain 
settings  of  the  correlator  made  for  the  smooth  surface  were  optimum  for 
the  rough  surface  case  as  well  and  no  further  adjustments  of  the  corre- 
lator or  oscilloscope  were  necessary  while  the  rough  surface  value  was 
being  computed.  In  a  few  cases  the  correlator  gain  was  changed  and  the 
ratio  of  the  two  gains  was  multiplied  times  the  computed  rough  surface 
correlation  before  the  relative  value  was  determined.  For  the  rough 
surface  case  the  amplitude  of  the  correlation  function  computed  was 
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proportional  to  the  mean  acoustic  pressure  received  at  the  hydrophone. 
C(t)  =  ^—   cos  («ot  -  *)   ,  (35) 


After  60  seconds  of  initial  integration,  the  peak-to-peak  value  of  this 
function  was  read  10  times  with  the  readings  taken  approximately  20 
seconds  apart.  The  average  value  of  these  readings  was  recorded  as  a 
"run".  The  reading  was  recorded  as  a  fraction  of  the  full  scale  reading 
and  represents  the  relative  mean  acoustic  pressure.  To  see  this  it  is 
noted  that  the  smooth  surface  correlation  function  was  set  to  equal  full 
scale  and  therefore  the  peak-to-peak  value  of  the  rough  surface  corre- 
lation function  is  equivalent  to  the  ratio: 

1<P^        -       |<P>|  (36) 


'20 


At  least  two  and  usually  three  runs  were  made  and  averaged  to  obtain 
each  data  point.  The  final  value  was  squared  to  yield  the  relative  mean 
acoustic  pressure  squared. 

A  procedure  was  followed  to  ensure  that  undesirable  electrical 
feed-through  of  the  driving  voltage  into  the  receiving  electronics  was 
not  introducing  error.  The  correlator  gains  and  oscilloscope  settings 
were  adjusted  to  get  optimum  correlator  performance  and  full  scale 
reading  for  the  smooth  surface  condition.  Then  without  making  further 
adjustments  the  aperture  was  blocked  so  that  no  sound  could  reach  the 
hydrophone.  The  cross-correlation  function  was  observed.  In  most  cases 
the  correlation  function  dropped  to  zero;  however,  in  the  parts  of  the 
experiment  where  the  low  source  level  Western  Electric  640AA  condenser 
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transmitter  was  used,  there  was  a  small  amount  of  coherent  electrical 
noise  detected  that  was  not  eliminated  from  the  system.  In  no  case  did 
this  exceed  .03  times  the  smooth  surface  correlation  amplitude, 
b.  Obtaining  the  Relative  Mean  Square  Acoustic  Pressure 

A  similar  procedure  was  used  to  compute  the  relative  mean 
square  acoustic  pressure.  In  this  case,  however,  an  auto-correlation  was 
performed  with  the  input  from  the  receiving  hydrophone  going  to  both 
correlator  channels.  Since  the  mean  square  value  of  a  signal  is  defined 
to  be  the  auto-covariance  evaluated  at  x  »  0  ,  the  output  of  the  corre- 
lator for  this  point  was  displayed  on  a  Varian  F-100  X-Y  Recorder. 
This  computation  does  not  eliminate  the  uncorrected  additive  noise 
introduced  by  the  hydrophone  and  the  receiving  electronics.  The  input 
can  be  expressed  as  e(t)  =  Po(t)  +  n(t)  and  auto-covariance  evaluated 
at  t  =  0  yields, 

R(0)  =  R2(0)  +  Rn(0)  (37) 


where  R.(0)  is  the  mean  square  acoustic  pressure  received  and  Rn(0) 
is  the  noise  power.  For  the  smooth  surface  condition, 

R20(0)  ■  p202  F2  D2  (38) 


In  the  experiment,  the  value  for  the  smooth  surface  mean 
square  acoustic  pressure  was  computed  first  and  plotted  against  time  on 
the  X-Y  plotter.  The  surface  was  roughened  and  the  correlator  output 
was  recorded  for  approximately  200  seconds  after  it  had  undergone  the 
initial  decay  from  the  smooth  surface  to  rough  surface  value  (approxi- 
mately 60  seconds).  The  average  value  during  the  "run"  was  determined 
visually  and  divided  by  the  smooth  surface  result. 
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R(0)    _-    R2(0)  +  Rn^Q)  (39) 

R0(0)       R2Q(0)  +  Rn(0) 


In  most  cases  in  the  acoustic  experiment  the  ratio 


Rn(0) 


was  less  than  .01  and  the  noise  terms  did  not  contribute  significantly 
to  the  relative  mean  square  values  obtained.  However,  in  some  cases, 
which  will  be  pointed  out  later,  tank  reverberation  became  a  source  of 
possible  error  which  could  not  be  ignored.  If  the  noise  terms  are 
neglected  then  the  expression  for  relative  mean  square  acoustic  pressure 
by  the  procedure  outlined  above  may  be  written  as 


<£,\         .  J^X        .  (40, 

\  '  P20  F  D2 
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IV.  PRESENTATION  AND  ANALYSIS  OF  THE  EXPERIMENTAL  DATA 

A.  THE  BASIC  ACOUSTIC  EXPERIMENT 

The  greatest  part  of  the  acoustic  data  was  obtained  using  the  basic 
acoustic  arrangement  and  the  two  well-studied  seas,  A  and  B.  The  source 
was  operated  in  four  positions,  either  directly  above  or  at  angles  up- 
wind of  the  aperture,  in  order  to  study  transmission  for  angles  of 
incidence  varying  from  normal  to  somewhat  greater  than  the  "critical 
angle"  for  plane  waves.  The  receiver  was  at  an  angle  either  in  the 
Snell  direction  or  near  grazing;  these  two  cases  were  believed  to  be  of 
particular  interest.  Five  frequencies  were  used,  ranging  from  10.0  kHz 
to  24.3  kHz.  In  terms  of  the  theoretical  roughness  parameter,  a  range 
of  roughness  from  R  =  .06  to  R  =  2.8  was  observed.  Study  of  higher 
roughness  was  precluded  by  the  limited  frequency  range  of  the  source 
used.  For  each  source-receiver  position,  (0-,  ,62)  ,  both  relative  mean 
coherent  intensity  and  total  intensity  (coherent  plus  incoherent)  were 
computed  with  the  correlator  for  each  sea  and  all  available  frequencies. 
The  resultant  data  have  been  plotted  against  the  roughness  parameter,  R, 
and  are  presented  in  Figures  12-17.  Figures  18-22  use  the  same  data 
plotted  against  R  but  each  graph  shows  the  relative  mean  acoustic  pres- 
sure squared  (coherent  intensity)  results  for  a  particular  frequency. 

1.  The  Relative  Mean  Coherent  Intensity 

Theory  predicts  that  the  coherent  component  of  the  received 

acoustic  intensity  will  decrease  exponentially  with  increasing  R  .  [R  = 

2  2  2 

k2  o  (c2/c,  cos  9-.  -  cos  e2)  ].  The  results  of  the  basic  experiment 

-R 
have  been  graphed  against  R  in  a  semilog  plot;  the  theoretical  e 

dependence  is  shown  on  each  graph.  The  first  observation  that  can  be 
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made  is  that  the  surface  appears  "rougher"  than  theory  admits.  Otherwise 
stated,  although  there  is  an  exponential  decrease  with  increasing  rough- 
ness, the  empirical  roughness  parameter  is  greater  than  the  theoretical 
value.  When  this  trend  was  first  noted  during  the  experiment  the 
reaction  was  to  reinstall  the  wave  probe  and  to  re-measure  the  mean 
square  wave  heights  of  the  two  seas  to  ensure  that  the  surface  had  not 
changed.  It  had  not. 

A  closer  examination  of  the  data  revealed  that  there  was  a 

distinct  difference  between  the  results  for  the  two  sea  surfaces.  The 

-ctR 
data  for  each  sea  appeared  to  fit  an  exponential  of  the  form  e    , 

where  a  is  a  constant  greater  than  unity  for  a  given  sea.  The  smaller 

sea,  A,  produced  a  greater  deviation  from  theory.  The  conclusion  was 

that,  unless  the  observed  results  were  products  of  experimental  errors, 

the  "acoustical  roughness"  of  these  surfaces  was  dependent  on  some  other 

parameter,  or  parameters  of  the  sea  surface  in  addition  to  the  mean 

square  wave  height.  However,  before  pursuing  this  interesting  hypothesis 

it  would  be  well  to  explore  the  sources  of  experimental  error. 

2.  Analysis  of  Experimental  Error 

There  were  a  number  of  factors  in  the  experiment  that  may  have 

introduced  uncertainty  in  the  results.  The  very  fact  that  the  process 

studied  is  a  stochastic  one  will  introduce  a  variance  dependent  on  how 

well  the  quantities  being  measured  have  been  averaged.  How  much  this 

variance  might  be  is,  of  course,  a  function  of  the  process  itself  and 

must  be  estimated  from  the  experimental  data.  As  noted  in  the  preceding 

section  on  procedures,  an  effort  was  made  to  ensure  that  the  technique 

used  was  giving  repeatable  results.  Since  each  run  consisted  of  a  series 

of  samples  of  the  values  computed  by  the  correlator,  a  sample  variance 
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could  be  calculated  and  a  confidence  interval  established.  But  the 
technique  used  in  this  processing  produced  an  unbiased  estimator  and, 
therefore,  this  source  of  uncertainty  seems  an  unlikely  candidate  to 
explain  a  consistent  trend  away  from  theory. 

The  uncertainty  in  the  determination  of  the  mean  square  wave 
height  might  appear  to  be  a  likely  candidate.  An  error  here  would  show 
up  as  a  trend  such  as  that  observed.  The  fact  that  the  two  surfaces,  and 
therefore,  the  two  determinations  of  mean  square  wave  height  gave  differ- 
ent results  reinforces  the  suspicion.  But  the  surface  was  measured  and 
re-measured  quite  carefully,  and  the  calibration  technique  seemed  ade- 
quate. Recalling  that  a  7.5%  error  in  the  mean  square  wave  height  was 
held  to  be  the  maximum  expected,  the  apparent  error  in  the  mean  square 

wave  height  can  be  calculated  and  compared.  For  sea  A  the  experimental 

-2R 
data  can  be  fit  approximately  to  e   .  This  would  indicate  a  mean 

square  wave  height  double  that  measured  or  a  100%  error.  Similarly  the 

apparent  error  for  sea  B  is  about  50%.  Clearly  the  observed  result  can 

not  be  explained  by  errors  in  wave  height  measurement. 

As  noted  in  the  discussion  of  the  anechoic  tank,  the  reverberant 
field  in  the  tank  was  not  insignificant.  An  analysis  of  the  possible 
effect  of  this  factor  is  in  order.  In  the  case  of  the  mean  acoustic 
pressure,  only  the  coherent  component  need  be  considered,  and  a  simpli- 
fied description  of  the  field  produced  at  the  receiver  can  be  used  to 
predict  the  effect  of  the  reverberant  field. 

We  consider  that  the  hypothetical  field  at  the  receiver  is  deter- 
mined by  the  vector  sum  of  components  from  three  sources:  1)  the 
ensonified  area,  2)  a  coherent  reflection  from  the  bottom  and  3)  a  second 
coherent  reflection  from  the  surface.  In  determining  the  relative  mean 
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field  the  smooth  surface  value  is  computed  first;  then  the  surface  is 
roughened.  The  contribution  direct  from  the  ensonified  area  is  reduced 
in  amplitude,  similarly  the  contribution  from  the  bottom  bounce  is 
reduced  in  amplitude,  although  perhaps  not  in  quite  the  same  proportion 
as  the  direct  contribution  because  of  the  geometry  factor  in  the  surface 
roughness  parameter.  But  a  more  significant  effect  may  be  caused  by  the 
third  component,  which  is  reduced  a  second  time  by  surface  scattering. 
This  should  be  most  noticeable  at  the  higher  frequencies,  where  the  loss 
due  to  surface  scattering  is  more  pronounced.  The  disproportionate 
reduction  in  the  amplitude  of  this  third  component  may  cause  a  larger 
or  smaller  result  for  the  relative  mean  field,  depending  on  the  relative 
phases  of  the  vectors  concerned. 

Although  the  preceding  argument  is  for  a  simple  case,  the  result 
should  hold  for  the  complicated  reverberant  field  that  actually  existed. 
This  effect  would  be  as  likely  to  produce  errors  in  one  direction  as  in 
the  other,  and  can  not  explain  a  consistent  trend  away  from  theory.  The 
net  effect  of  reverberation  must  be  to  introduce  some  uncertainty  into 
the  experimental  results.  That  this  has  caused  some  problems  is  evident, 
particularly  for  the  higher  frequencies  and  near  horizontal  transmission. 
For  example,  Figure  22  gives  the  results  for  24.3  kHz  and  the  inconsis- 
tencies of  these  data  are  quite  apparent.  For  the  (6-.  ,02)  combinations, 
(12°, 80°),  (16°, 80°),  and  (16°, 85°)  the  roughness  parameter,  R,  is  about 
the  same.  But  there  was  large  variation  observed  in  the  relative  co- 
herent intensity.  This  effect  can  be  noted  for  both  the  "high"  and  the 
"low"  sea  but  is  more  pronounced  for  the  "high"  sea.  It  is  not  sur- 
prising that  reverberation  becomes  a  more  significant  factor  for  this 
particular  geometry.  Consider  the  field  plot  for  24.3  kHz  and  8-j  =  16° 
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field  the  smooth  surface  value  is  computed  first;  then  the  surface  is 
roughened.  The  contribution  direct  from  the  ensonified  area  is  reduced 
in  amplitude,  similarly  the  contribution  from  the  bottom  bounce  is 
reduced  in  amplitude,  although  perhaps  not  in  quite  the  same  proportion 
as  the  direct  contribution  because  of  the  geometry  factor  in  the  surface 
roughness  parameter.  But  a  more  significant  effect  may  be  caused  by  the 
third  component,  which  is  reduced  a  second  time  by  surface  scattering. 
This  should  be  most  noticeable  at  the  higher  frequencies,  where  the  loss 
due  to  surface  scattering  is  more  pronounced.  The  disproportionate 
reduction  in  the  amplitude  of  this  third  component  may  cause  a  larger 
or  smaller  result  for  the  relative  mean  field,  depending  on  the  relative 
phases  of  the  vectors  concerned. 

Although  the  preceding  argument  is  for  a  simple  case,  the  result 
should  hold  for  the  complicated  reverberant  field  that  actually  existed. 
This  effect  would  be  as  likely  to  produce  errors  in  one  direction  as  in 
the  other,  and  can  not  explain  a  consistent  trend  away  from  theory.  The 
net  effect  of  reverberation  must  be  to  introduce  some  uncertainty  into 
the  experimental  results.  That  this  has  caused  some  problems  is  evident, 
particularly  for  the  higher  frequencies  and  near  horizontal  transmission. 
For  example,  Figure  22  gives  the  results  for  24.3  kHz  and  the  inconsis- 
tencies of  these  data  are  quite  apparent.  For  the  (0,  ,0p)  combinations, 
(12°, 80°),  (16°, 80°),  and  (16°, 85°)  the  roughness  parameter,  R,  is  about 
the  same.  But  there  was  large  variation  observed  in  the  relative  co- 
herent intensity.  This  effect  can  be  noted  for  both  the  "high"  and  the 
"low"  sea  but  is  more  pronounced  for  the  "high"  sea.  It  is  not  sur- 
prising that  reverberation  becomes  a  more  significant  factor  for  this 
particular  geometry.  Consider  the  field  plot  for  24.3  kHz  and  8-j  =  16° 
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of  Figure  8  ;  the  SPL  at  85°  is  6  dB  less  than  that  at  30°  which,  when 
scattered  from  the  bottom,  no  doubt  produced  a  proportionately  large 
reverberant  field  at  the  receiver.  It  seems  reasonable  to  conclude  that 
if  the  data  appear  well  behaved  and  consistent  then  the  effects  of  the 
reverberant  field  must  be  small;  but  when  the  data  become  inconsistent, 
with  large  swings  in  the  experimental  value  for  small  changes  in  rough- 
ness, the  reverberation  must  be  introducing  significant  errors.  Despite 
these  apparent  errors  in  some  parts  of  the  experiment,  the  fact  remains 
that,  in  general,  there  was  an  unexplained  trend  that  clearly  indicated 
that  the  experimental  acoustical  roughness  is  greater  than  theory 
predicts. 

3.  The  Relative  Mean  Total  Intensity 

Theory  predicts  that  for  low  roughness  the  incoherent  component 
of  the  acoustic  intensity  is  dominated  by  the  coherent  component  and, 

therefore,  the  relative  mean  square  acoustic  pressure  should  follow  an 

-R 

e   dependence.  As  R  increases  the  mean  square  value  should  depart 

-R 
from  e   and,  for  transmission  in  the  Snell  direction  at  very  high 

roughness,  the  relative  mean  acoustic  intensity  should  be  proportional 

to  R~  .  In  this  experiment  R  has  been  limited  to  values  less  than  3.0 

and  reverberation  effects  further  restricted  the  range  of  roughness 

where  reliable  results  have  been  obtained.  In  general,  it  would  seem 

prudent  to  dismiss  as  unreliable  the  mean  square  data  for  those  cases 

which  showed  significant  reverberation  effects  in  the  relative  coherent 

intensity.  When  the  more  reliable  data  is  considered  it  can  be  seen 

that  the  qualitative  predictions  of  theory  are  correct,  at  least  over 

the  limited  range  of  R  that  has  been  observed. 
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B.  THE  SUPPLEMENTARY  ACOUSTIC  EXPERIMENTS 

The  validity  of  the  theory  within  the  limitations  imposed  by  the 
model  (and  therefore  for  realistic  "real  world"  conditions)  must  be 
seriously  questioned  based  on  the  results  of  the  basic  experiment.  It 
is  true  that  the  qualitative  behavior  of  rough  surface  transmission 
appeared  to  fit  the  description  of  the  theory  but  it  is  also  true  that 
the  surfaces  were  considerably  "rougher"  acoustically  than  predicted. 
There  were  a  number  of  assumptions  and  approximations  that  were  required 
to  derive  the  theoretical  results  and  the  model  experiment  was  designed 
with  the  knowledge  that  some  of  these  assumptions  were  incorrect  and 
some  of  the  approximations  were  questionable. 

The  following  are  the  significant  violations  of  theoretical  con- 
straints: 

1)  An  anisotropic  sea  with  a  non-Gaussian  spatial  correlation 
function  and  a  near  Gaussian  but  skewed  wave  height  probability 
distribution  was  used. 

2)  The  incident  wave  was  not  plane. 

3)  The  requirements  of  the  Kirchhoff  approximation,  L  »  ^  ,  were 
not  satisfied. 

4)  For  the  incoherent  component,  ^fA  was  not  much  greater  than  L. 
Of  these,  the  third  violation  was  considered  to  be  the  most  serious. 
The  first  seemed  less  critical  because  the  only  specification  of  the 
rough  sea  in  the  derivation  of  the  mean  acoustic  field  requires  that  the 
distribution  of  heights  be  Gaussian  and  this  was  nearly  so.  The  non- 
plane  wave  front  was  not  considered  to  be  the  primary  cause  for  the 
observed  deviation  but  it  was  desired  to  eliminate  this  possible  cause 
in  order  to  isolate  other  effects.  The  fact  that  different  seas  gave 


85 


different  results  in  the  basic  experiment  motivated  the  attempt  to  per- 
form some  supplementary  experiments  that  would  investigate  the  effect  of 
changing  the  sea  surface  "shape". 

1.  Dependence  on  Sea  Surface  Characteristics 

The  first  supplementary  experiments  to  be  tried  were  with  the 
newly  generated  and  measured  seas  C  and  E.  All  other  components  of  the 
basic  experiment  were  unchanged  but  as  described  in  Section  III,  Sea  C 
had  a  greater  correlation  length  than  A  but  had  approximately  the  same 
mean  square  wave  height.  Sea  E  was  a  very  small  sea  (short  correlation 
length)  that  violated  the  Kirchhoff  assumption  more  strongly  than  did 
sea  A. 

In  all  the  supplementary  experiments  only  the  normal  incidence 
case  with  the  hydrophone  in  the  Snell  direction  (e?  =  0°)  was  studied. 
In  the  basic  experiment  normal  incidence  appeared  to  be  the  least 
affected  by  reverberation. 

Figure  23  compares  the  results  obtained  for  seas  A  and  C.  Sea 
C  appears  "smoother"  acoustically.  If  the  upwind  downwind  correlation 
distance,  L,  is  taken  to  be  representative  of  the  minimum  radius  of 
curvature  of  the  surface,  then  the  Kirchhoff  approximation  assumes  that 
LA?  »  1  .  For  sea  A  the  ratio  LuAo  is  less  than  or  equal  to  .35 
and  for  sea  C,  LUA2  <   .52  .  The  obvious  conclusion  is  that  sea  C  has 
produced  results  closer  to  theory  because  the  "shape"  of  the  surface  was 
more  nearly  like  the  class  of  surfaces  for  which  the  theory  was  derived, 
that  is  L  A?  is  larger  at  each  frequency  (or  roughness). 

The  results  obtained  with  sea  E,  shown  in  Figure  24,  were  ob- 
tained by  a  slightly  different  technique  than  that  outlined  in  the  data 
collection  procedures  section.  Since  very  small  transmission  losses  were 
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FIGURE  23.  RELATIVE  COHERENT  INTENSITY  VS  R 
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FIGURE  24.   RELATIVE  TOTAL  INTENSITY  VS  R 
SEA  E 
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involved,  the  relative  mean  square  acoustic  pressure  was  assumed  to  be 
equal  to  the  relative  mean  acoustic  pressure  squared  and  the  latter  value 
was  estimated  by  computation  of  the  autocovariance  at  t  =  0  .  This 
technique  permitted  greater  resolution  since  the  quantity  calculated  was 
the  square  of  a  value  slightly  less  than  one,  and,  therefore,  is  easier 
to  differentiate  from  one  than  the  value  itself.  In  addition  the  com- 
puted value  could  be  plotted  on  an  X-Y  recorder  and  read  with  more 
accuracy. 

In  spite  of  this  approach,  experimental  uncertainties  must  have 
been  relatively  large  when  compared  to  the  small  transmission  losses. 
With  these  reservations  noted  it  can  be  observed  that  the  deviation  from 
theory  does  appear  to  be  somewhat  greater  than  that  seen  for  sea  A.  For 
sea  E  the  maximum  L  /Ap  ratio  is  .13  .  Another  way  of  expressing  the 
results  obtained  with  sea  E  is  that  for  the  very  small  seas,  where  the 
Kirchhoff  approximation  is  most  strongly  violated,  the  transmission  loss 
is  so  small  that  the  absolute  error  introduced  by  the  violation  of  theory 

doesn't  amount  to  much.  However,  for  sea  E  which  had  the  greatest  devi- 

-2  4R 
ation  from  the  Kirchhoff  approximation,  the  empirical  slope,  e  *   , 

is  furthest  from  theory. 

Since  the  constraint  of  the  Kirchhoff  approximation  is  expressed 

in  terms  of  the  ratio,  LA2  ,  it  was  desirable  to  extend  the  experiments 

to  higher  frequencies  with  the  same  seas.  This  would  supplement  the 

data  of  the  basic  experiment  by  increasing  the  range  of  roughness  observed 

and,  it  was  hoped,  would  provide  some  more  information  on  the  observed 

deviation  from  theory.  In  order  to  do  this  without  the  influence  of  the 

non-plane  wave,  a  parabolic  reflector  source  was  used.  This  permitted 

extension  of  the  experiment  up  to  43.0  kHz,  while  maintaining  a  nearly 
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plane  incident  wave.  The  source  was  also  used  with  a  new  sea  that  was 
generated  with  the  longest  possible  correlation  length  and  a  moderate 
mean  square  wave  height. 

2.  Dependence  on  Wavefront  Curvature 

The  plane  wave  results  (obtained  with  the  parabolic  reflector 
source)  are  given  in  Figures  25,  26  and  27  where  they  are  compared  with 
data  from  Figure  12.  An  additional  source  of  experimental  error  was 
encountered  in  the  plane  wave  experiments.  A  small  but  not  insignificant 
electrical  feedthrough  of  correlated  noise  could  not  be  eliminated.  The 
relative  amplitude  of  this  interference  was  estimated  as  described  in 
Section  III-B,  and  the  error  flags  in  the  plots  correspond  to  the  cases 
where  this  component  either  added  in  phase  or  180°  out  of  phase.  The 
data  points  without  error  flags  were  obtained  with  no  detectable  elec- 
trical feedthrough. 

Comparing  the  results  obtained  with  seas  A  and  B  and  the  plane- 
wave  with  those  for  the  diverging  wave  (Figures  25,  26),  it  appears  that 
there  is  no  appreciable  difference  due  to  wave  front  curvature.  This 
would  seem  to  indicate  that  wavefront  curvatures  at  least  as  large  as 
.17A-,  across  the  ensonified  area  do  not  significantly  effect  the  coherent 
pressure  received  and  that  a  Fresnel  correction  is  not  required. 

3.  Dependence  on  Sea  Surface  Correlation  Length 

The  great  advantage  of  a  model  experiment  is  that  it  can  move 
into  areas  relatively  unexplored  by  theorists  and,  by  studying  the  effects 
of  various  parameters  in  the  model  gain  insights  about  a  process  that  may 
be  quite  difficult  to  describe  analytically.  The  parameter  that  seemed 
to  have  a  significant  effect  not  predicted  by  theory  was  the  correlation 
length,  L.  The  correlation  length  was  measured  for  all  sea  surfaces  used, 
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FIGURE  25.  RELATIVE  COHERENT  INTENSITY  YS  R 
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O  Plane  Wave 
+  Diverging  Wave 


0.1 


0.5 


1.0 


1.5 


91 


0.1 


0.01 


FIGURE  26.  RELATIVE  COHERENT  INTENSITY  VS  R 
SEA  B   e1  =  0°  ,  e„  =  Oc 
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FIGURE  27.  RELATIVE  COHERENT  INTENSITY  VS  R 
SEA  D   0-,  =  0°  >  e2  =  0° 
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and  there  seems  to  be  a  relationship  between  the  correlation  length  and 
the  amount  of  deviation  from  theory. 

Although  the  data  that  can  be  used  to  make  an  empirical  deter- 
mination of  this  apparent  relationship  are  limited,  there  seems  to  be 
enough  to  derive  an  approximate  correction  factor  to  be  applied  to  the 
theoretical  predictions.  The  correction  factor  developed  here  uses  the 
approximation  that  each  sea  with  correlation  length  L  and  L  /x?  <  1  , 
will  produce  a  relative  mean  acoustic  pressure  squared  that  can  be  pre- 
dicted by 

where  a  =  a(L  )  .  The  model  study  did  not  extend  beyond  the  regime 
Lu/A2  <  1  . 

In  order  to  determine  an  empirical  function  for  a(L  )  ,  the 
normal  incidence  data  from  the  model  study  was  used  to  estimate  a  value 
for  a  for  each  sea  used.  A  best  fit  straight  line  on  a  semilog  plot  was 
drawn  for  each  sea  as  shown  in  Figures  23-27.  Assuming  the  theory  would 
hold  for  large  correlation  lengths  required  that  a(L  )  approach  one  as 
L  became  large.  This  led  to  the  speculation  that  a(L  )  could  be 
expressed  as  follows: 


a(Lu)  =  1  +  be"Lu/Lo 

where  b  and  L  are  constants.  To  determine  the  values  for  these 
constants  and  to  see  if  the  experimental  data  actually  did  fit  such  a 
function,  the  experimental  values  of  a  -  1  were  plotted  against  the 
experimental  correlation  length  and  compared  to  the  function: 


a  -  1  =  be 


Lu/Lo 
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Since  seas  A,  B  and  D  were  the  most  thoroughly  studied,  particular 
attention  was  made  to  fit  the  empirical   relationship  to  them.     A  good 
fit  was  obtained  as  shown  in  Figure  28.     The  empirical   correction  to 
theory  thus  obtained  is  given  by 

a0-u)     =     1+4.8  e"Lu/1,75  [Lu  in  cm] 

In  order  to  demonstrate  the  effect  of  this  added  correction 
factor  all  the  normal  incidence  relative  coherent  intensity  data  have 
been  plotted  against  aR  ,  the  modified  roughness,  in  a  form  normalized 
with  the  value  predicted  by  the  corrected  theory.  It  can  be  seen  from 
Figure  29  that  the  bulk  of  the  data  plots  are  within  10%  of  the  pre- 
dicted value.  This  scatter  is  easily  attributed  to  the  experimental 
uncertainties. 

It  can  be  noted  that  for  large  values  of  aR  Figure  29  shows 
a  rapid  departure  from  the  corrected  theory.  This  phenomenon  is  reminis- 
cent of  an  effect  reported  by  Mayo  [8]  and  starts  at  about  the  same  value 
of  roughness.  Figures  26  and  27  also  show  this  trend  as  a  leveling  off 
of  the  transmission  loss  of  the  coherent  intensity  at  the  higher  fre- 
quencies. This  might  have  been  interpreted  as  a  trend  toward  theory  as 
the  L  A2  ratio  becomes  larger.  It  is  noted,  however,  that  sea  B 
(Figure  26)  departs  from  the  empirical  roughness  prediction  between 
R  =  2.5  and  R  =  3.0  ,  where  L  A?  =  .635  ,  while  sea  D  (Figure  27) 
departs  between  R  =  3.0  and  R  =  3.5  ,  where  LuA2  =  1.2  .  It  would 
appear  that  there  is  little  connection  between  the  observed  "leveling 
off"  of  the  coherent  intensity  and  the  L  A«  ratio.  The  significant 
similarity  between  both  of  these  cases  is  that  they  appear  to  depart 
from  theory  at  approximately  the  same  value  of  aR,  '■>-'4.0  . 
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FIGURE  28.   EMPIRICAL  DETERMINATION  OF  a(L 

(Points  corresponding  to  Seas  A-E 
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The  corrected  theory  can  be  used  to  predict  the  behavior  of  the 
incoherent  intensity  as  well  as  the  coherent  intensity.  The  relative 
incoherent  intensity  can  be  obtained  by  subtracting  the  relative  coherent 
intensity  (33)  from  the  relative  total  intensity  (34).  In  the  present 
experiment,  unfortunately,  the  difference  of  two  nearly-equal  and  some- 
what uncertain  quantities  must  be  taken.  Specializing  the  theory  for 
transmission  in  the  Snell  direction,  introducing  the  correction  factor, 
and  writing  in  a  convenient  form  yields: 


<^A 


nC0^erent      =  S(aR)  (41) 

Tfl//A 

where  S(aR)  is  the  function  given  in  (28). 

Recall  that  the  correlation  length,  L,  referred  to  in  the  theory 
is  that  for  an  isotropic  surface  with  a  Gaussian  correlation  function. 
The  factor  ttL  /A  is,  in  effect,  a  ratio  of  "correlated  area"  to  ensoni' 
fied  area.  To  compare  the  results  obtained  with  the  anisotropic  sea 
surfaces  one  uses  L  =  L  L   where  L  and  L  are  the  upwind  and  cross- 
wind  correlation  lengths,  respectively.  Using  the  Snell  angle  data  for 
seas  A  and  B  and  (e^e^  angles  (0°,0°),  (6°, 30°),  (12°, 60°),  the  ratio 
(41)  is  plotted  against  aR  in  Figure  30. 

Comparison  with  the  predicted  S(aR)  behavior  shows  reasonable 
agreement  over  the  range  of  aR  observed.  Since  the  value  of  S(aR)  is 
changing  slowly  with  aR  over  the  region  (aR  =  0.5  to  5)  where  the 
incoherent  component  was  determined,  equally  good  agreement  would  have 
been  obtained  if  the  data  were  plotted  against  R  and  compared  to  S(R). 
The  higher  roughnesses  will  have  to  be  studied  to  obtain  the  complete 
picture  of  behavior  of  the  incoherent  component. 
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4.  A  Point  Source  Experiment 

A  simplified  model  of  the  ocean-acoustic  experiments  conducted 
with  the  FLIP  was  set  up  in  the  laboratory.  The  objective  of  this  model 
experiment  was  to  gain  some  insights  that  might  aid  in  the  analysis  of 
the  at-sea  experiment. 

The  source  used  was  the  Western  Electric  640AA  transmitter 
located  one  meter  above  the  unmasked  water  surface.  A  detailed  discus- 
sion of  the  source  and  acoustic  arrangement  is  given  in  Section  III-B. 
The  hydrophone  arrangement  was  the  same  as  that  in  the  basic  experiment. 
It  was  located  one  meter  from  the  point  on  the  surface  directly  below 
the  source  and  could  be  positioned  accurately  at  any  angle  from  this 
point.  The  filter  of  the  Hewlett-Packard  Model  3590A  Wave  Analyzer  was 
used  to  provide  a  10  Hz  bandwidth  for  the  smooth  surface  data  and  a  100  Hz 
bandwidth  for  the  mean  square  acoustic  pressure  measurements. 

Medwin  [11]  has  developed  an  expression  for  the  underwater  sound 
field  of  a  point  source  in  air  above  a  smooth  surface.  The  result  was 
derived  using  a  ray  theory  approach  and  gives  the  received  acoustic  pres- 
sure amplitude,  p?.  ,  in  terms  of  the  pressure  amplitude,  p-,0  ,  on  the 
surface  directly  beneath  the  point  source.  The  geometry  of  the  problem 

can  be  seen  in  the  following  figure. 

source 
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The  hydrophone  position  was  measured  in  terms  of  92Q  ,  while  the  pre- 
diction is  expressed  in  terms  of  e2  as  well  as  e, ,  H,  and  D. 
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The  point  source  model  was  used  to  experimentally  determine 
smooth  surface  underwater  acoustic  field.  The  field  was  plotted  by 
sweeping  e2Q  (and  the  hydrophone)  from  0°  to  180°  while  the  hydrophone 
output  was  filtered  and  detected  by  the  wave  analyzer  and  then  plotted 
in  dB  vs.  62Q  .  A  simple  iterative  calculation  with  the  IBM  360/67 
Computer  was  used  to  determine  9?  from  e2Q  so  that  the  experimental  data 
could  be  easily  compared  with  the  prediction.  Reasonable  agreement  was 
obtained  for  the  relative  sound  pressure  levels  for  10.0  and  20.0  kHz  as 
a  function  of  6«q  as  shown  in  Figures  31  and  32.  The  absolute  SPL's 
are  not  compared  in  these  figures  because  of  the  possible  interference 
of  the  reverberant  field. 

A  second  probe  of  the  underwater  sound  field  as  a  function  of 
depth  for  a  frequency  of  20  kHz  produced  a  plot  of  SPL  relative  to  that 
at  the  surface  that  can  be  compared  to  the  prediction  of  Medwin's  formula 
with  more  confidence.  This  plot  is  given  in  Figure  33,  where  the  solid 
line  indicates  the  predicted  field.  The  increased  variation  of  the 
experimental  SPL  with  depth  is  attributed  to  the  reverberant  field.  The 
experimental  and  predicted  fields  are  in  substantial  agreement  and  appear 
to  be  approaching  a  spherical  spreading  dependence  with  depth  as  would 
be  expected  in  the  "far  field". 
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Some  data  were  taken  for  the  rough  surface  (sea  D)  sound  trans- 
mission from  a  point  source.  For  this  experiment  the  hydrophone  was 
positioned  at  62q  =  0°  and  a  depth  of  one  meter.  The  relative  mean 
total  intensity  was  measured  in  the  manner  described  previously  and  the 
results  are  shown  in  Figure  34.  Interpretation  of  these  results  is  not 
easily  done;  a  quantitative  prediction  to  compare  with  these  data  is 
not  available  at  this  time.  Figure  34  shows  the  relative  mean  square 
pressure  beginning  to  fall  with  increasing  roughness  but  then  leveling 
off  with  a  transmission  loss  of  less  than  1  dB.  Qualitatively  this 
effect  can  be  attributed  to  incoherent  contributions  to  the  field  at 
the  receiver  from  areas  that  do  not  contribute  significantly  to  the 
smooth  surface  field.  A  quantitative  analysis  of  this  phenomenon  awaits 
further  theoretical  work  as  well  as  more  extensive  model  studies. 
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V.  CONCLUSIONS 

The  model  sea  experiment  has  revealed  that  the  theory  based  on  the 
Kirchhoff  method  of  evaluating  the  Helmholtz  integral  is  not  correct 
for  many  of  the  surfaces  and  acoustic  frequencies  that  are  most  important 
in  the  problem  of  air-ocean  sound  transmission.  However,  the  theory, 

which  predicts  that  the  coherent  component  of  the  transmitted  intensity 

-R 
will  decrease  exponentially  (e  )  with  the  roughness  parameter,  R,  can 

be  used  if  an  empirically-derived  correction  factor  is  introduced.  The 

correction  factor  is  a  function  of  the  surface  correlation  length  only 

and  multiplies  the  theoretical  roughness  parameter,  R,  to  give  a  corrected 

acoustical  roughness,  aR.  Defining  the  "correlation  length", Lu,  as  the 

minimum  distance  over  which  the  sea  surface  spatial  correlation  falls 

to  e~  ,  the  correction  factor,  a,  reduces  to  unity  for  L  A„  >>  1  and 

corrects  for  the  case  where  the  sea  surface  correlation  length  is  less 

than  or  equal  to  the  acoustic  wavelength  in  water. 

When  the  correction  factor  is  applied,  the  experimental  results  agree 
with  the  prediction  of  coherent  sound  transmission  over  a  range  from  low 
to  moderate  acoustic  roughness  but  depart  significantly  from  the  corrected 
theory  for  the  roughness  greater  than  aR  -  4.0  .  This  departure  may  be 
related  to  a  similar  phenomenon  that  has  been  observed  for  the  coherent 
component  of  rough  surface  scattering  by  Mayo  [8]  and  others. 

The  predictions  of  the  corrected  theory  seem  to  describe  the  behavior 
of  the  total  acoustic  intensity.  The  incoherent  component  of  the 
received  signal  appears  to  follow  the  predicted  behavior  within  the 
limits  of  experimental  accuracy.  Experimental  limitations  prevented  the 
testing  of  theory  for  the  ^ery   high  roughness  case. 
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The  wavefront  curvature  does  not  appear  to  be  a  critical  factor  for 
the  range  of  curvature  studied.  Comparison  of  the  results  of  the  same 
experiment  made  with  a  near-plane  wave  generated  by  a  parabolic  reflector 
and  a  diverging  wavefront  with  as  much  as  .17x,  curvature  across  the 
ensonified  area  showed  no  detectable  difference. 

The  limited  experiments  conducted  with  a  nearly  omnidirectional 
point  source  above  an  unlimited  water  surface  seem  to  verify  a  theoretical 
prediction  for  the  smooth  surface  underwater  sound  field  derived  by 
Medwin.  The  rough  surface  transmission  results  could  be  qualitatively 
explained  based  on  theory. 

A  quantitative  prediction  for  the  point  source  problem  may  be 
achieved  by  partitioning  the  ensonified  area  into  subareas  that  meet 
the  specifications  of  the  limited  ensonified  area  used  in  the  theory  and 
by  summing  the  contributions  of  each  subarea.  This  procedure  may  not 
prove  too  difficult  for  determining  the  coherent  component  but  it  is 
not  a  trivial  extension  of  the  theory  in  the  apparently  important  case 
of  the  incoherent  contributions  to  the  received  intensity;  for  the 
incoherent  contributions  from  subareas,  the  surface  height  correlation 
between  subareas  must  be  accounted  for. 

The  discovery  that  correlation  length  must  be  used  to  empirically 
correct  the  roughness  parameter  seems  to  be  consistent  with  the  fact  that 
the  experiment  strongly  violated  the  Kirchhoff  assumption,  L/x^  >>  1. 
It  suggests  that  additional  theoretical  work  needs  to  be  done  for  the 
case  of  surfaces  and  acoustic  frequencies  that  do  not  meet  the  Kirchhoff 
restriction.  The  question  arises  as  to  whether  the  size  of  the  ensoni- 
fied area  with  respect  to  the  correlation  length  might  become  a  factor 
for  the  coherent  component  as  it  is  in  the  present  theory  for  the 
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incoherent  component  of  intensity.  The  answer  to  this  and  other  questions 
brought  up  by  this  study  must  be  found  before  this  wery   interesting 
problem  can  be  considered  to  be  solved. 
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APPENDIX  A  -  STUDY  OF  SEA  SURFACE  CHARACTERISTICS 

1.  Probability  Density  Functions  of  Surface  Wave  Heights 

The  probability  density  of  wave  heights  was  determined  for  the  wind- 
driven  seas  A  and  B  as  part  of   a  study  which  used  the  CI5000/XDS9300 
hybrid  computer  to  process  twenty-minute  wave  records  of  each  sea.  The 
procedures  used  in  this  study  were  as  follows: 

a.  The  wave  probe  was  located  at  the  center  of  the  ensonified  area 
and  oscillated  to  provide  a  10-minute  calibration  signal  to  the  FM  tape 
recorder.  After  the  calibration  signal,  twenty-minute  wave  records  of 
each  sea  were  recorded. 

b.  The  calibration  signal  and  the  wave  records  were  replayed  at 
ten  times  speed  into  the  analog-to-digital  converter  of  the  hybrid 
computer  facility.  The  converter  sampled  the  wave  record  at  an  effective 
real-time  rate  of  50  times  per  second  until  an  ensemble  of  25,600  values 
was  obtained. 

c.  The  ensemble  mean  was  computed  and  then  adjusted  to  zero  by 
subtracting  the  computed  mean  from  each  sample  value.  The  mean  square 
voltage  for  the  two  wave  records  and  the  calibration  signal  was 
computed  and  the  calibration  was  used  to  determine  root  mean  square  wave 
heights,  0  • 

d.  The  digital  computer  then  sorted  the  ensemble  values  into 
intervals  of  width  0.2a  .  The  results  were  normalized  to  give  unit 
area  and  are  plotted  in  Figures  Al  and  A2  .  .  A  Gaussian  distribution 
curve  is  included  for  comparison. 

2.  Spatial  Correlation  of  the  Sea  Surface 

The  two-dimensional  correlation  function  of  the  sea  surfaces  were 
determined  by  cross-correlation  of  the  wave  record  from  the  two  wave 
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probes,  A  and  B,  and  recording  the  value  computed  with  zero  time  delay. 
The  wave  probes  were  mounted  on  traverse  that  permitted  one  probe  to  be 
positioned  at  the  center  of  the  ensonified  area  and  the  other  moved  in 
any  radial  direction  to  any  separation  within  about  30  cm.  The  wave 
probes  were  not  calibrated  as  the  correlation  function  could  be  computed 
as  follows 


CU)  =  


[CAA(T-"0)  WT=0)] 


1/2 


where  £  is  the  probe  separation. 

In  order  to  determine  the  correlation  for  a  given  separation  it  was 
necessary  to  record  the  cross-correlation  on  the  x-y  recorder  for  at 
least  200  seconds  and  visually  average  the  resulting  plot.  The  mean 
square  voltage  of  each  probe  was  computed  with  a  probe  separation  of 
6  cm.  and  used  to  normalize  the  cross-correlations  obtained  for  separations 
greater  than  4  cm.  For  separations  of  4  cm.  or  less  it  was  necessary 
to  determine  the  mean  square  voltage  for  each  probe  at  each  separation. 
This  was  required  because  the  probe  sensitivity  decreased  with  separation 
until  a  spacing  of  at  least  3  to  4  cm.  was  reached.  Because  of  the  long 
averaging  times  required  to  compute  each  point  of  the  Spatial  Corre- 
lation Function  this  was  an  extremely  time  consuming  process. 

Spatial  correlation  was  determined  for  the  upwind  and  crosswind 
directions  as  well  as  four  intermediate  angles.  Figures  A3  -  A7  show 
the  correlation  functions  obtained. 
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3.  Wave-Number  Sea  Surface  Spectrum 

The  two-dimensional  surface  spatial  correlation  data  was  used  to 
compute  a  wave-number  spectrum  for  sea  A.  The  two-dimensional  fourier 
transform  of  the  correlation  function,  C(£,n)  yields  the  normalized 
directional  spectrum  of  the  sea  surface: 

-   B(t)  -  (  ^  )2   //  C(5.n)  e-i(t?)  d? 

where    r  =  i,£  +  i'2n         ,       K  =  i,K    +  i'2K      and     |fc|   =  2tt/A       , 

A  =  surface  wavelength. 

To  accomplish  this  transformation  the  HARM  subroutine  of  the  IBM/360 

Scientific  Subroutine  Package  was  used.     This  subroutine  uses  the  Cooley- 

Tukey  (1965)   fast  fourier  transform  algorithm.     The  spatial  correlation 

data  was  formed  into  a  16  by  16  array  of  sample  values,   (£,n)»  with  one 

centimeter  spacing.     The  upwind  direction  was  taken  to  be  the  x-dimension 

of  the  array.     The  transformation  of  this  array  by  HARM  yielded  an  8  by  8 

array  of  fourier  coefficients  for  wave  numbers   (K   ,K  ). 

x  y 

The  resolution  of  the  spectrum  was  limited  by  the  maximum  separation 
obtained  in  the  surface  spatial  correlation.  Although  it  would  have  been 
desirable  to  extend  the  correlation  function  out  to  greater  distances, 
the  uncertainty  in  such  data  due  to  the  limited  averaging  time  of  corre- 
lator would  have  made  the  accuracy  of  the  resultant  spectrum  somewhat 
questionable.  If  the  correlation  function  had  been  sampled  every  0.5  cm 
the  spectrum  could  be  extended  to  higher  wave-numbers.  This  appears 
justified;  however,  it  was  not  done. 

The  advantage  of  a  wave-number  spectrum  as  a  description  of  a  sea 
surface  is  that  it  clearly  demonstrates  the  directional  nature  of  the  sea. 
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Figures  A8,  A9,  and  A10  show  the  wave-number  spectrum  in  the  upwind- 
downwind  direction  (K  =0),  the  cross-wind  spectrum  (K  =0)  and  45°  off 

y  x 

the  wind   (K  =K  ). 
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APPENDIX  B  -  SUMMARY  TABLES  OF  DATA 
1.  Basic  Experiment  Data 
COHERENT  AND  TOTAL  INTENSITY: 


SEA  A  (a  =  .173  cm,  Lu  =  2.1  cm) 


freq 

el 

e2 

10.0 

0 

0 

13.5 

0 

0 

17.5 

0 

0 

20.3 

0 

0 

24.3 

0 

0 

10.0 

6 

30 

13.5 

6 

30 

17.5 

6 

30 

20.3 

6 

30 

24.3 

6 

30 

10.0 

12 

60 

13.5 

12 

60 

17.5 

12 

60 

20.3 

12 

60 

24.3 

12 

60 

10.0 

12 

80 

13.5 

12 

80 

17.5 

12 

80 

20.3 

12 

80 

24.3 

12 

80 

10.0 

16 

80 

13.5 

16 

80 

17.5 

16 

80 

20.3 

16 

80 

24.3 

16 

80 

10.0 

16 

85 

13.5 

16 

85 

17.5 

16 

85 

20.3 

16 

85 

24.3 

16 

85 

<p>l  <pp*> 


.059  .89  .90 

.107  .77  .78 

.180  .69  .68 

.240  .57  .60 

.342  .52  .54 

.063  .81  .82 

.115  .72  .78 

.193  .62  .64 

.259  .53  .55 

.366  .48  .50 

.074  .85  .82 

.135  .71  .72 

.23  .57  .58 

.306  .56  .60 

.43  .40  .40 

.087  .90  .90 

.16  .67  .75 

.27  .64  .69 

.35  .55 

.51  .36  .49 

.085  .84  .81 

.154  .72  .68 

.26  .67  .64 

.345  .45 

.49  .29 

.088  .88  .94 

.16  .74  .79 

.27  .56  .62 

.36  .67 

.51  .41  .42 
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COHERENT  AND  TOTAL  INTENSITY  (Continued) 
SEA  B  (a  =  .41  cm,  l_u  =  3.5  cm) 


freq 

ei 

e2 

R 

l<p>l2 

<PP*> 

10.0 

0 

0 

.33 

.60 

.63 

13.5 

0 

0 

.60 

.39 

.44 

17.5 

0 

0 

1.01 

.20 

.25 

20.3 

0 

0 

1.35 

.15 

.18 

24.3 

0 

0 

1.91 

.078 

.091 

10.0 

6 

30 

.353 

.50 

.53 

13.5 

6 

30 

.635 

.34 

.38 

17.5 

6 

30 

1.12 

.18 

.26 

20.3 

6 

30 

1.45 

.12 

.19 

24.3 

6 

30 

2.04 

.060 

.15 

10.0 

12 

60 

.42 

.52 

.57 

13.5 

12 

60 

.75 

.32 

.37 

17.5 

12 

60 

1.32 

.15 

.20 

20.3 

12 

60 

1.71 

.12 

.20 

24.3 

12 

60        * 

2.4 

.056 

.18 

10.0 

12 

80 

.49 

.64 

.64 

13.5 

12 

80 

.88 

.28 

.40 

17.5 

12 

80 

1.50 

.17 

.25 

20.3 

12 

80 

2.0 

.11 

-- 

24.3 

12 

80 

2.82 

.038 

.12 

10.0 

16 

80 

.475 

.57 

.59 

13.5 

16 

80 

.86 

.31 

.38 

17.5 

16 

80 

1.45 

.192 

.24 

20.3 

16 

80 

1.95 

.144 

-- 

24.3 

16 

80 

2.75 

.022 

-- 

10.0 

16 

85 

.50 

.60 

.68 

13.5 

16 

85 

.89 

.31 

.42 

17.5 

16 

85 

1.51 

.16 

.24 

20.3 

16 

85 

2.02 

.19 

-- 

24.3 

16 

85 

2.86 

.09 

.19 
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INCOHERENT  INTENSITY: 

freq 

61 

e2 

R 

aR 

<pp*>incoh. 

\pp  /incoh. 
ttL2/A 

SEA  A 

(a  =   2.1    , 

,  ttL2/A 

=    .06) 

10.0 

0 

0 

.059 

.123 

.01 

.17 

13.5 

0 

0 

.107 

.224 

.01 

.17 

17.5 

0 

0 

.180 

.378 

.00 

.00 

20.3 

0 

0 

.240 

.50 

.03 

.50 

24.3 

0 

0 

.342 

.72 

.02 

.34 

10.0 

6 

30 

.063 

.132 

.01 

.17 

13.5 

6 

30 

.115 

.232 

.06 

1.0 

17.5 

6 

30 

.193 

.405 

.02 

.336 

20.3 

6 

30 

.259 

.543 

.02 

.336 

24.3 

6 

30 

.366 

.77 

.02 

.336 

10.0 

12 

60 

.074 

.150 

.00 

.00 

13.5 

12 

60 

.135 

.284 

.01 

.17 

17.5 

12 

60 

.23 

.483 

.01 

.17 

20.3 

12 

60 

.306 

.642 

.04 

.67 

24.3 

12 

60 

.43 

.90 

.00 

.00 

SEA  B 

(a  =  1.43, 

,  nt2/A 

=    .21) 

10.0 

0 

0 

.33 

.47 

.03 

.14 

13.5 

0 

0 

.60 

.86 

.05 

.24 

17.5 

0 

0 

1.01 

1.44 

.05 

.24 

20.3 

0 

0 

1.35 

1.93 

.03 

.14 

24.3 

0 

0 

1.91 

2.72 

.013 

.054 

10.0 

6 

30 

.353 

.505 

.03 

.14 

13.5 

6 

30 

.635 

.405 

.04 

.19 

17.5 

6 

30 

1.12 

1.59 

.08 

.38 

20.3 

6 

30 

1.45 

2.08 

.07 

.33 

24.3 

6 

30 

2.04 

2.92 

.09 

.43 

10.0 

12 

60 

.42 

.60 

.05 

.24 

13.5 

12 

60 

.75 

1.07 

.05 

.24 

17.5 

12 

60 

1.32 

1.88 

.05 

.24 

20.3 

12 

60 

1.71 

2.45 

.08 

.38 

24.3 

12 

60 

2.4 

3.43 

.124 

.59 
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2.  Supplementary  Experiment  Data 


freq      R       aR 
SUPPLEMENTARY  SEAS: 


■aR 


I(P> 


-aR 


SEA  C     (a  =    .17  cm,   Lu  =  3.2  cm,  a  =   1.52)      (6^0°,  62=0< 


10.0 

.059 

.09 

.914 

.88 

.962 

13.5 

.107 

.163 

.850 

.82 

.965 

17.5 

.18 

.274 

.760 

.78 

1.03 

20.3 

.24 

.365 

.696 

.71 

1.02 

24.3 

.342 

.52 

.595 

.58 

.975 

SEA  E 

(a  =  .051 

cm,  Lu  = 

.8  cm,  a  = 

3.7) 

(6^0°, 

e2=o°) 

10.0 

.0049 

.018 

.982 

.987 

1.01 

13.5 

.0089 

.033 

.968 

.973 

1.01 

17.5 

.0151 

.056 

.946 

.935 

.99 

20.3 

.0203 

.075- 

.928 

.955 

1.03 

24.3 

.029 

.107 

.899 

.94 

1.04 

PLANE  WAVE   (PARABOLIC  REFLECTOR  SOURCE): 


SEA  A 

(a  =  .173 

cm,  Lu 

=  2.1  cm,  a  =  2.1) 

(0^0°, 

e2=o°) 

10.0 

.059 

.124 

.884 

.93 

1.05 

13.5 

.107 

.225 

.798 

.80 

1.00 

17.5 

.18 

.378 

.686 

.71 

1.03 

20.3 

.24 

.505 

.604 

.68 

1.13 

22.3 

.293 

.615 

.541 

.48 

.89 

24.3 

.342 

.72 

.488 

.51 

1.04 

27.0 

.43 

.90 

.406 

.42 

1.03 

30.0 

.53 

1.11 

.33 

.29 

.88 

34.5 

.70 

1.47 

.230 

.225 

.98 

37.0 

.81 

1.70 

.183 

.176 

.96 

43.0 

1.09 

2.29 

.109 

.109 

1.00 

SEA  B 

(a  =  .41 

cm,  Lu  = 

3.5  cm,  a  =  1.42) 

(6^0°, 

e2=o°) 

10.0 

.33 

.468 

.626 

.63 

1.01 

13.5 

.60 

.852 

.427 

.39 

.915 

17.5 

1.01 

1.43 

.240 

.24 

1.00 

20.3 

1.35 

1.92 

.147 

.16 

1.09 

22.3 

1.63 

2.32 

.098 

.096 

.98 

24.3 

1.91 

2.72 

.066 

.069 

1.05 

27.0 

2.40 

3.41 

.033 

.038 

1.15 

30.0 

2.96 

4.20 

.015 

.0256 

1.70 

34.5 

3.91 

5.55 

.0039 

.0132 

3.4 

37.0 

4.5 

6.4 

.0017 

.008 

4.7 
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l(p)l 

e-R 


PLANE  WAVE   (PARABOLIC  REFLECTOR  SOURCE)      (Continued) 
freq  R  aR  e"aR  |<p)|2 

SEA  D     (a  =   .32  cm,  Ly  =  4.2  cm,  a  =   1.26)      (0^0°,   62=0°) 

10.0  .20  .252  .778  .81  1.04 

13.5  .365  .46  .632  .562  .89 

17.5  .617  .777  .465  .482  1.03 


2Q.0 

.80 

1.01 

.366 

.372 

1.02 

24.0 

1.15 

1.45 

.235 

.230 

.98 

27.0 

1.46 

1.84 

.160 

.162 

1.01 

30.0 

1.80 

2.27 

.104 

.090 

.87 

34.5 

2.38 

3.0 

.050 

.050 

1.00 

37.0 

2.74 

3.45 

.032 

.029 

.91 

43.0 

3.69 

4.65 

.0096 

.015 

1.56 

LARGER  APERTURE   (50  x  50  CM): 

freq  R     '  |(p)|2  <pp*> 

SEA  A     (a  =  0.173  cm,  Lu  =  2.1   cm)     (0^0°  ,62=0°) 


10.0 

.059 

.79 

.78 

13.5 

.107 

.61 

.58 

17.5 

.18 

.62 

.67 

20.3 

.24 

.59 

Not  done 

24.3 

.342 

.41 

.47 

SEA  B  (a  •■ 

=  0.41  cm,  Lu  = 

3.5  cm) 

(e, 

=o°,e2=o°) 

10.0 

.331 

.405 

.44 

13.5 

.595 

.23 

.34 

17.5 

1.01 

.23 

.63 

20.3 

1.35 

.15 

Not  done 

24.3 

1.91 

.04 

.425 

(Note:     The  larger  aperture  data  shows  significant  interference 
from  reverberation  —  it  was  not  used  in  the  analysis.) 
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3.  Point  Source  Experiment  Data 
SEA  D  (a  =  .32  cm,  L„  =  4.2  cm) 


freq 

R 

<PP*> 

(S/N)1 

(S/N)2 

10.0 
13.5 
17.5 
20.0 
24.0 

.20 
.365 
.617 
.80 
1.15 

.95 
.81 
.85 
.86 
.87 

20  dB 
20  dB 
15  dB 
15  dB 
14  dB 

30  dB 
30  dB 
25  dB 
25  dB 
25  dB 

Note: 

(S/N)1 

is  signal 

to 

total 

noise  ratio 

(S/N)2  is  correlated  signal  to  noise  ratio  and  is  a  measure 
of  the  coherent  electrical  feedthrough. 
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The  corrected  theory  appears  to  be  valid  for  both  coherent  and  incoherent 
components  of  intensity  for  all  values  of  L/A2  provided  that   R  <  4.0  .  For 
higher  values  of  roughness  there  is  a  "leveling  off"  of  the  transmission  loss 
of  the  coherent  component  similar  to  that  recently  observed  but  unexplained  in 
rough  surface  scattering. 


DD 


FORM      l/fTQ 

S/N    0101 -807-681  1 


(PAGE     1 ) 


133 


Security  Classification 


Security  Classification 


Sound  Transmission 
Underwater  Sound 
Air-Sea  Interface 
Rough  Water  Surface 
Wave  Height  Distribution 
Wave  Spectrum 
Correlation 
Coherent 


DD,F„r..1473  'back, 


S/N  0101 -807-682* 


134 


Security  Classification 


liord  


SHEtF   BINDER 

^SS     Syrocuse,  N.  Y. 
— —     Stockton,  Colif. 


I 


J  SEP  71 


2  5  5  3  2 


The 

H13J 
c 


122456 


Thesis  ._?7I 

H1135         Hagy 

c.l         Transmission  of 
sound  through  a 
randomly  rough  air- 
sea  interface. 

•vi  /fas  c/ 

s  1 
2  i 

2553  2 


.5  JAN  60 


Thesis 

H1135 

c.l 


Hagy 

Transmission  of 
sound  through  a 
randomly  rough  air- 
sea  interface. 


122456 


thesH1135 

Transmission  of  sound  through  a  randomly 

llllliilill 

3  2768  001  03720  3 

DUDLEY  KNOX  LIBRARY 


^^ 


